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ABSTRACT 


The computer program AFTBDY generates a body-fitted 
curvilinear coordinate system for a wedge curved after body. 
(Figure 1). This wedge curved after body is being used in an 
experimental program at the Langley Research Center. The 
coordinate system generated by AFTBDY will be used to solve 
3D compressible N.S. equations. The coordinate system in 
the physical plane is a cartesian x,y,z system, whereas, in 
the transformed plane a rectangular ?,rifC system is used. 

The coordinate system generated is such that in the transformed 
plane coordinate spacing in the direction is constant 

and equal to unity. The physical plane coordinate lines in 
the different regions are clustered heavily or sparsely depending 
on the regions where physical quantities to be solved for by the 
N.S. equations have high or low gradients. The coordinate dis- 
tribution in the physical plane is such that x stays constant in 
n and 5 direction, whereas, z stays constant in ^ and n direction. 
The desired distribution in x and z is input to the program. 
Consequently, only the y-coordinate is solved for by the program 


AFTBDY. 



List of Symbols 


D 

D* 

IMAX 

JiyiAX 

KMAX 

J 

P,Q,R 

S.O.R 

x,y,z 


e 



0 ) 


Subscripts 
?/ri/C,x,y,z 
CC ,nn ,xx,yy etc. 
5Ti/xy,xz etc. 

Superscripts 

(s) 


Physical Plane 
Transformed Plane 

Maximum number of points in ? - direction 

Maximum number of points in n - direction 

Maximum number of points in 5 - direction 

Jacobian 

Inhomogeneous terms used for coordinate attraction 

Successive-Over-Relaxation iteration 

Physical coordinates 

Transformed coordinates 

Error criterion for convergence 

Boundary contour surfaces in the physical plane 


Boundary contour surfaces in the transformed plane 


Acceleration parameter for Gauss-Seidel iteration 


Denotes first partial differentiation 
Denotes second partial differentiation 
Denotes corss partial differentiation 

Denotes iteration nximber 



I. Numerical Generation of Body-Fitted 
Curvilinear Coordinate Systems 

The method for numerically generating the Boundary-Fitted 
curvilinear coordinate system is presented as applicable in 
the case of the wedge-curved after body. The coordinate system 
generated is to be used to solve 3-D compressible Navier Stokes 
equations. The coordinate system in the physical plane is a 
cartesian x,y,z system, whereas in the transformed plane rectangu- 
lar system is used. The coordinate system generated is 

such that in the transformed plane, coordinate spacing in all 
three C/Hf? direction is constant and equal to unity. The- physical 
plane coordinate lines in the different directions are clustered 
heavily or sparsely depending on the regions where the physical 
quantities to be solved by the N.S. equations have large or small 
gradients. The coordinate distribution in the physical plane is 
such that X stays constant in n and C direction, whereas z stays 
constant in 5 and n direction. The desired distribution in x and 
z direction is input to the program. Consequently, only the y 
coordinate is solved for by the program AFTBDY. The clustering 
of points in y direction is specified by a point distribution on 
C = 1, ? = 1 line. 

Section A provides a brief discussion of the mathematical 
formulation and section B describes the numerical technique used 
to generate the coordinate system. 



A. Mathematics of Transformation 


Consider the transformation of a simply connected 3-d 

region D (Figure 2), into a right angle parallepoid region D* 

(Figure 3). We require that faces T, , F., F_, F- map 

****** ±^J45o 

. F F F F F F ... 

into i' 2'’3'4'5'6 ^©spectively . For identification 

purposes, region D is referred to as the physical plane, and 
* 

D as the transformed plane. In the D plane it is assumed 
that X is constant for a specified n and ? and z is constant 
for a specified C and n. In the present problem only y is 
solved for since the desired distribution in x and z direction 
is input to the program. A 3-d formulation from x,y,z to C,D/? 
coordinates is defined as follows: 
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specified functions. 


The inhomogeneous terms P,Q, and R are selected to control 
the spacing of n = constant and C = constant lines in the physical 
plane and several forms for these terms can be used. In the pre- 
sent case a special form of P, Q and R was used based on the point 
distribution on ab, ac, and ae lines (Figure 2) . (For convenience 
Co = n-, = C, =1 and Cc “ IMAX, ri- = JMAX and ?. = KMAX is chosen 
thus ensuring a field size of IMAX x JMAX x KMAX) . 

In order that the transformed plane contain a uniform grid, 
the dependent and independent variables must be interchanged in 
equations (2a) , (2b) and (2c) . The resulting equations are given by, 
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The functions a, b, c, d, e, £ are specified by the known 

surfaces F. , F., F., F_ and F-. On all six surfaces F, / 

b 1 ^ 

r^, F^, Fg and Fg Dirchlet boundary conditions are used. 

As indicated before/ in the present problem only y is solved 
for since the desired distribution in x and z is input to the pro- 
gram. The last aspect to be discussed is hov/ P, Q and R are computed. 
P is computed based on x variation of points on line ab (Figure 2) . 
Similarly, Q and R are computed based on y, z variation on line ac 
and ae, respectively. The following discussion will make this point 
clearer. 

If one assumes on ab that the y and z variation is zero (y^ = 
z^ = 0) , on ac x and z variation is zero (x^ = z^ = O) , and on ae x 
and y variation is zero (x = y^ = 0) then from equations (4) and 
(5) we obtain 



_ _ 
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Notice P, Q, 

and R are functions of C,n, 

and Zf respectively. 



Thus only one dimensional arrays of size IMAX, JMAX, KMAX are re- 
quired to store P, Q, and R, respectively. 
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B. Numerical Solution Techniques 


Based on the above discussion equation (4b) must be 
numerically solved. The coefficients are given by 

equations (5) and P, Q, R are given by equations (7) . The 
numerical solution of the transformation can be easily per- 
formed. First P, Q, and R are computed as follows: 

^ x(I+l>l,l) - 2x(I,l,l) + x(I-l,l>l) 

(x(I+l>l,l) - x(I-l,l,l) )**3 (Sa) 

2^I<_IMAX-1 

Q(j) = y(l.J+l,l) - 2y(l,J>l) + y(l,J-l.l) 

(y(l,J+l,l) - y(l,J-l,l) )**3 

2<^J^f^lC-l (8b) 

R(K) = z(l,l,K+l) - 2z(l,l,K) + z(l,l,K-l) 

(z(l,l,K+l) - z(l,l,K-l) )**3 

2<_K<_MAX-1 (8c) 

A S.O.R Gauss Seidell iteration is used to numerically 
solve equation (4b) . All derivatives in (4b) are written in 
second order central difference form. Then following expression 
for intermediate value of y at point (I, J, K) is obtained, 

y(I,J,K) = CF*[c^^(y(I+l,J,K) + y(I-l,J,K)) - 2 

+2Ci3y^^ + C22(y(l/J+1^K) + y(I,J-l,K)) 
+C 33 (y(I,J,K+l) + y(I,J,K-D) - 2 C 23 y^^ + 
J^(P(I)y^ + Q(J)y^ + R(K)yj,)] (Sa) 
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where 

Cp = 0.5/(c^j^ + C 22 + C 23 ) (9b) 

In (9a) we have 2<_I^IMAX-1, and 2^K^KMAX-1. 

In addition Cj ^2 etc. are given by equation (5) . If f 

denotes x, y, or z then the derivatives of k, Yt z in difference 


form required for equation (9a) are given by, 

f^ = 0.5(f (I+1,J,K) - f(I-l,J,K)) (10a) 

f^ = 0.5(f (I,J+1,K) - f(I,J-l,K)) (10b) 

f^ = 0.5(f (I,J,K+1) - f(I,J,K-D) (10c) 

f^^ = 0.25(f (I+1,J+1,K) - f(I-l,J+l,K) 

-f (I+1,J-1,K)+ f (I-1,J-1,K) ) (lOd) 

f = 0.25(f (I,J+1,K+1) - f(I,J-l,K+l) 

n ^ 

-f (I,J+1,K-1) + f (I,J-1,K-1) ) (lOe) 

f^^ = 0.25(f (I+1,J,K+1) - f(I-l,J,K+l) 

-f (I+1,J,K-1) + f (I-1,J,K-1) ) (lOf) 

f^^ = f(I+l,J,K) - 2f(I,J,K) + f(I-l,J,K) (lOg) 

f^^ = f(I,J+l,K) - 2f(I,J,K) + f(I,J-l,K) (lOh) 

f^^ = f(I,J,K+l) - 2f(I,J,K) + f(I,J,K-l) (lOj) 


If (9a) is the intermediate value for S.O.R iteration then 


value of y at the next iteration (S + 1) is as follows: 



( 11 ) 


(I,J,K) = (1-w) y^®^(I,J,K) + w y(I,J,K) 

In equations (9) and (10) most current iterative value of y are 
used and u is the acceleration parameter. 

In equations (9) and (10) most current iterative value of y are 
used and w is the acceleration parameter. 

The convergence is reached whenever the absolute value of 
y(I, J, K)^® y(I, J, K) is less than a specified quantity 

£. Generally a value of 0.000005 to 0.00001 is specified for e. 

A value of 0.4 to 0.8 is used for u. 

C. Initial Guess and Boundary Point Specification 

The input to the program consists of x coordinate specifica- 
tion on line ab, y specification on line ac, and z specification 
on line ae. This implies that {x(I,l,l), l£l£lI1AX} , {y(l,J,l), 
l£J_<JMAX} and {z(l,l,K), l£K£KMAX} are known. (The program is 
written for a body for which x varies from o to 90.0, and z varies 
from 0 to 12.0). The initial guess is specified as follows, 

x(I,J,K) = x(I,l,l) (12a) 

z(I,J,K) = z(l,l,K) (12b) 

where l^I^IMAX, 1^J<JI4AX, l£K<KMAX. 

The y coordinates are specified as follows. On surface Fg 
(J=JMAX) the coordinates are given by, 

y(I,JMAX,K) = y(l,JMAX,l) C13a) 


KKIMAX and KK<KMAX. 



(13b) 


on r 


also 


and 


y on 


2 ( J =1) y coordinate is specified as follows: 

y(I,l,l) = x(I,l,l) * tan (4°) 

if x(I,l,l) < 62.3 

y(I,l,l) = C+[D^ - (x(I,l,l) - E)^] **0.5 

if x(I,l,l)> 62.3 

C = -457.8350 
D = 463.3164 

E = 94.56704 

y(I,l,K) = y(I,l,l) 

if x^ 62 

/ 

y(I,l,K) = y(I,l,l) - R + [r2 - z^j 
R = [A + B(x(I,l,l) “ 62.3)]**0.5 
A = 0.00035521 
B = 0.000223505 
2^I<IMAX and 2<_K^KMAX. 

remaining of the points is given as follows; 

y(I,J,K) = y(I,l,K) + (y (I, ■ 

(y(l,J,l) - y(l,l,l))/ 
y(l,JMAX,l) - y(l,l,D) 
1<KIMAX, 2<J<JMAX-1 and 1<K<KI4AX 


(13c) 

3 (13d) 

if x>62.3 (13e) 

y(I,l/K))* 

(13f) 



line of symmetry 





a* e* f* b* is surface aefb in figure 2. 


longitudinal curvature 


spanwise curvature 


Figure 1. Wedge Curve After Body 
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surface bfhd 
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Physical Plane 
Figure 2 . (continued) 
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surface a'e'g'c' 
surface a'e'f'b* 
surface e'f'h'g' 
surface b'f'h'd' 
surface c'g'h'd' 


Transformed Plane 


Figure 2. Coordinate Transformation 



II. Computer Program 


A. Program Listing 

A computer program to numerically generate boundary-fitted 
coordinate system for a wedge after body combination was written 
in FORTRAN IV language to run on the LaRC CDC NOS system complex. 

The computer program consists of the main program AFTBDY and sub- 
routines INIT, CALCOR, IPRTC, and IPLTC. A brief description of 
the main program and the four subroutines is provided in Chapter IV. 
The program listing follows this page. 



CREATION RUN CAROS ENCOUNTERED IN INPUT ’ UPOATE 1.2-452 

♦DECK AFTBOr 

PROGRAM AFTBOY(INPUT#OUTPUTfSOLNXX»TAPEl-SnLHXX#TAPE5»INPUT» 

1 TAPE6-0UTPUTI 


C^^ ♦ 

C*^ THIS PROGRAM GENERATES COORDINATE SYSTEM FOR A 3-0 AFTER ♦ 

C*+ BODY. THE COORDINATES GENERATED ARE BASED ON THE METHOD ♦ 

€♦♦ DEVELOPED BY THOMPSON# THAMES AND HASTIN OF MISSISSIPPI ♦ 

C+* STATE UNIVERSITY, ♦ 

C + % ♦ 

C++ FOR FURTHES ENQUIRIES CONTACT i “ ♦ 

C + + ♦ 

C*+ DR, DILIP KUMAR • ♦ 

C++ DR. JULIUS HARRIS ♦ 

C++ NASA - LANGLEY RESEARCH CENTER ♦ 

C++ HAIL STOP 163 ’ / * ♦ 

C + + HAMPTON# VA 23665 ' '♦ 

C++ PHONEI 804-827-3696 ♦ 

C++ FTS CODE 928-1110 " ♦ 

C + + . ♦ 


c ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦+♦♦+♦♦♦ ♦♦♦♦♦♦♦♦♦♦♦%♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

C + + 

C++ THIS PROGRAM GENERATES COORDINATES FOR A AFTER BODY. DUE 

C++ TO SYMMETRY ONLY HALF OF THE BODY IS CONSIDERED. THE 

C++ COORDINATES IN X AND 7 DIRECTION ARE FIXED. SOLUTION IS 

C++ OBTAINED FOR Y ONLY. THE DESIRED DISTRIBUTION OF POINTS 

C++ IS DEFINED ON J-1#K-1 AND I ■ 1#IMAX POINTS# THEN ON 

C++ I«1#K-1 AND J»1#JMAX AND FINALLY ON I»1#J«1 AND K«1#KMAX. 

C++ THIS DISTRIBUTION CAN BE SPECIFIED BASED ON UNIFORM# 

C++ EXPONENTIAL# OR POLYNOMIAL OIS TR IPUTION. BASED ON THIS 

C++ DISTRIBUTION THE ATTRACTION COEFFICIENTS P#Q# AND R ARE 

C++ COMPUTED, THE DISTRIBUTION OF POINTS MUST BE SPECIFIED 

C++ SO THAT P#Q# AND R ARE CQNTTNOUS, BASED ON THESE P#Q#AND R 

C++ THEN Y COORDINATE IS COMPUTED, 

C + + 

C + + CARD INPUT 

C + + 

C++ CAROl t SPECIFIES GRID SIZE 

C++ +♦ IHAX# JMAX#KMAX (FORMAT 315) 

C + + IMAX I NO OF NODES IN X-DIRECTION 

C++ JHAX I NO OF NODES IN Y-DIRECTION 

C + + KMAX I NO OF NODES IN Z-DIRECTION ' 

C + + 

C + + CARD2 : GAUSS-SEIOELL P AR , #CONVEPGENCE CRITERION# ERROR 
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EATION RUN 

C + + 
C + * 

C* + 
€♦♦ 

C + + 
C + + 
C++ 
C + + 
C + + 
C++ 
C + + 
C++ 
C++ 
C + + 
C + + 
C+ + 
C++ 


CAROS ENCOUNTERED IN INPUT ' “ UPDATE K2-A52;~ 

PRINT FREQUENCY FLAG# AND MAXIMUM NO OF ITERATIONS 
♦♦ R1#R2#R3#R^#IERR#ITERHX (FORMAT 4F10*0#2I5) 

R1 I GAUSS SEIDELL ITERATION PARAMETER 

(USUALLY 0.6) 

R2#R3#i MAXIMUM ERROR ALLOWED IN X#Y#AND I DIRECTION ' ' 

RA RESPECTIVELY (USUALLY 0.000005 TO 0.00001) 

lERR : INTERVAL AT WHICH MAX ERROR IS PRINTED. 

ITERMXl MAXIMUM NUMBER OF ITERATIONS PERFORMED. 


CARD3 t 
++ 


CAR04 t 
♦ + 


BOUNDARY POINT COORDINATES 
X#YfZ (FORMAT 3F10.0) 

THERE ARE IMAX+JMAX+KMAX SUCH CAROS. IMAX CARDS ARE 
INPUT FOR AND !•) 10 IMAX. SIMILARY CAROS 

FOR I«l>K-liJ»l TO JMAX AND I*liJ«l»K«l TO KMAX ARE 
INPUT. FOR IMAX CAROS ONLY X IS INPUT# 

FOR JMAX CARO Y# AND KMAX CAROS ONLY 1 IS INPUT. 

OTHER COORDINATES ARE LEFT BLANK AND ARE COMPUTED 
BY THE PROGRAM. 

PRINT OPTION PARAMETERS 

IPRTl#IDELlf JDEL1#IPRT2#!0EL2# JDEL2 (FORMAT 615) 
IPRTl I PRINT INITIAL GUESS O-NO 1-YES. 

IDELl I INTERVAL OF I AT WHICH INITIAL GUESS IS 

PRINTED. 1 MEANS FOR EVERY 1# 2 MEANS FOR !• 
1#3#5#7 AND SO ON. 

JDELl : INTERVAL OF J-COORDINATE FOR WHICH INITIAL 
GUESS IS PRINTED. 

1PRT2 I SAME AS I PRT 1 # IDELl # JOELl FOR CONVERGED OR 
IDEL2 PARTIALLY CONVERTED SOLUTION. PARTIALLY 
JDEL2 CONVERGED SOLUTION IS PRINTED IN CASE THE 
SOLUTION DOES NOT CONVERGE AFTER ITERMX 
ITERATIONS. 


CARDS t HEADER CARO 

HED1(I)#I«1#A (FORMAT ^AlO) 

++ HEDl t 40 BYTES OF HEADER TNFO PRINTED ON COORDINATE 
PRINTOUT AND PLOTS 


CAR06 t ADDITIONAL HEADER CARD 

HED2(D# I-l#4 ( FORM AT4 AID) 

SAME AS CARDS. (BOTH CARDS AND CARD6 REQUIRED) 


CARD? I INITIAL GUESS PLOT OPTION 
NPLTl (FORMAT 15) 


79/06/22 


AFTBDY 

AFTBOY 

AFTBDY 

AFTBDY 

AFTBOY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 
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AFTBDY 

AFTBOY 

AFTBDY 
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CREATION RUN 


CAROS ENCOUNTERED IN INPUT 


update 1,2-A52, 


79/06/22. 18. 29,34 


CAA 

♦♦ NPLTl » SPECIFTES NUMBER 

OF SURFACES FOR WHICH 

AFTBDY 

89 

CAA 

INITIAL GUESS IS 

PLOTTED. 

AFTBDY 

90 

CAA 



AFTBDY 

91 

CA* 

CARDS t there are NPLTl CARD TYPE 8. THESE CARDS SPECIFY 

AFTBDY 

92 

CA* 

VARIOUS SURFACES THAT NEED TO BE PLOTTED. 

AFTBDY 

93 

CA# 

ITYPE#ISTN#SIZEl 

(FORMAT RISiFlO.OJ 

AFTBDY 

94 

CAA 

♦♦ ITYPE 1 TYPE OF SURFACE 

.-11 STATION. - 2 K STATION 

AFTBDY 

95 

CA* 

ISTN : VALUE OF T OR K 

STATION DEPENDING ON ITYPE. ‘ 

AFTBDY 

96 

CA* 

SIZEl 1 MAX PLOT SIZE IN 

Z DIRECTION IF ITYPE ■ 1# 

AFTBDY 

97 

CA* 

OR MAX PLOT SIZE 

IN X-DIRECTION IF ITYPE • 2'. 

AFTBDY 

98 

CA + 



AFTBDY 

99 

CA* 

CARD9 t SIMILAR TO CARO TYPE 7 

AMO 8 RESPECTIVELY FOR 

AFTBDY 

100 

C*A 

CARDIO CONVERGED OR PARTIALLY CONVERGED SOLUTION. 

AFTBDY 

101 

CA* 



AFTBDY 

102 

CAA 



AFTBDY 

103 

CAA 

AAA SAMPLE RUNSTREAM 


AFTBDY 

104 

CAA 



AFTBDY 

105 ' 

CAA 



AFTBDY 

106 

CAA 

AFTB0Y#T7777#CM320000. 

BIN 34 D KUMAR 

AFTBDY 

107 

CAA 

USER/479019C. 


AFTBDY 

108 

CAA 

CHARGE>100718fLRC. 


AFTBDY 

109 

CAA 

GcTCQLDPL»AFTBDY/UN«375732NJ. 


AFTBDY 

110 

CAA 

UPDATE(F>C«VISFL1). 


AFTBDY 

111 

CAA 

COPYBP# INPUT#EOTFTl. 


AFTBDY 

112 

CAA 

REWIND#VISFL1#£DTFT1. 


AFTBDY 

113 

CAA 

EDIT# VISFLlfl-EDTFTl. 


AFTBDY 

114 

CAA 

REWIND#V1SFL1. 


AFTBDY 

115 

CAA 

FTN(PL»77777>I-VISFL1). 


AFTBDY 

116 

CAA 

at.tach>lrcgosf/un«library»na. 


AFTBDY 

” 117 

CAA 

LOS ET^LIB-LRCGOSF* PRES ET-ZERO. 

‘ 

AFTBDY 

118 

CAA 

LGO. 


AFTBDY 

— 119 

CAA 

PLOT.VARIAN 


AFTBDY 

120 

CAA 

REW1ND#S0LNXX. 


AFTBDY 

121 

CAA 

C0PY8R#S0LNXX#AFS0LN. 


AFTBDY 

122 

CAA 

REWINO#AFSQLN. 


AFTBDY 

123 

CAA 

SAVE# AFSOLN. 


AFTBDY 

124 

CAA 

7/8/9 


AFTBDY 

125 

CAA 

AIDENT DUMl 


AFTBDY 

126 

CAA 

AINSERT AFTBDY.4 


AFTBDY 

127 

CAA 

CAA 


AFTBDY 

128 

CAA 

7/8/9 


AFTBDY 

129 

CAA 

RSW MXLL/#/ 26/;a 


AFTBDY 

130 

CAA 

RSi/ MXLL/#/ 26/|A 


AFTBDY 

131 

CAA 

RSI/ MXLL/#/ 21/lA 


AFTBDY ' 

132 



RUN 

€♦♦ 

C^* 

C*+ 

C++ 
C++ 
C++ 
C + + 
C + + 
C + + 
C + + 
C+ + 
C + + 
C + + 
C + + 
C++ 
C + + 
C++ 
C + + 
C + + 
C + + 


C+ + 


CARDS ENCOUNTERED IN INPUT ' UPDATE 1.2-A52; — 


END 

7/0/9 

INPUT DATA, SEE EXPLANATION ABOVE, 

6/7/8/9 

NOTEl AFTER FIRST 7/8/9 CARD 3 CARDS ARE INPUT TO 
SATISFY INPUT FOR UPDATE UTILITY, 

AFTER 2ND 7/B/9 CARO EDIT CARDS ARE INPUT, THE 
PROGRAM IS WRITTEN SO THAT MAXIMUM NUMBER OF “ 

POINTS IN X,Y*Z DIRECTION CAN BE CHANGED BY SIMPLE 

INPUT. IN THE ABOVE CASE FIELD SIZE OF 26#26#21 

IS ASSUMED. IN THE THREE CARDS AFTER RSi/ AND MXLL 
I#J#K SHOULD BE INPUT RESPECTIVELY. THIS COULD NOT 
BE SHOWN ABOVE BECAUSE TEXT EDITOR CHANGES TEXT IN 
COMMENTS CAROS ALSO, 

IT IS ASSUMED THAT PROGRAM RESIDES ON FILE AFTBDY 
UNDER USER NO. OF 375732N AND USER ^79019C HAS 
PERMISSION TO USE IT, FURTHERMORE# IT IS ALSO ASSUMED 
THAT THE COORDINATE SYSTEM SOLUTION IS TO BE STORED 
ON FILE AFSOLN, THIS NAME CAN BE CHANGED BY 
CHANGING THREE CAROS TOWARDS THE END OF THE RUNSTREAM. 

THE PLOTS ON CALCOMP CAN BE OBTAINED BY REPLACING 
PLOT.VARIAN CARO BY FOLLOWING CARDS, 
PL0T.CALP0ST#ll(X0«1.0# Y0»0,2 ) 

CONT. //PAPER DO# 

CONT, LEROY PEN SIZE 0.3 

CONT. MULTIPLE PLOT MODE 

CONT. THESE PLOTS APE FOR REPORTS// 


COMMON /BXXX/ X ( IMXLL# JMXLL# KMXLU # Y ( IMXLL# JMX LL# KMXLL ) # 

1 Z(IMXLL#JMXLL#KMXLL)#P( IMXLL) #0(JMXLL)»P (KMXLU 
COMMON /AXXX/ IMAX# JMAX#KMAX# ILIMIT# JL IMTT#KL IMIT# ISQLN# 

1 KKK#XERRMX#YERRMX#7EPRMX,IXERI»IXERJ#IXERK#IYERI# lYERJ# 

2 lYERK# IZERI# IZERJfl ZERK,Rl#R2#P3#RA#IERR# ITERMX#IPRT1# 

3 IDEL1#JDEL1# IPRT2# I 0EL2# J0EL2# NPLTl # ITYPE# ISTN,SIZE1# 

4 SIZE2#LL#HED1C4)#HE02C4)#IPLSET#FILLXXC25) 

DATA KKK# IXER I# I XER J# I XERK # I YER I # I YER J# I YERK# I ZER I » I ZER J# IZERK# 
1 XERRMX# YERRMX#ZERRMX#IPLSET/10+0#3+0,0#0/ 

ILIMIT ■ IMXLL 
JLIMIT • JMXLL 


79/06/22, 

18,29 

AFTBDY 

133 

AFTBDY 

134 

AFTBDY 

135 

AFTBDY 

136 

AFTBDY 

137 

AFTBDY 

13B 

AFTBDY 

139 

AFTBDY 

140 

AFTBDY 

141 

AFTBDY 

142 

AFTBDY 

143 

AFTBDY 

144 

AFTBDY 

145 

AFTBDY 

146 

AFTBDY 

147 

AFTBDY 

148 

AFTBDY 

149 

AFTBDY 

150 

AFTBDY 

151 

AFTBDY 

. 152 

AFTBDY 

153 

AFTBDY 

154 

AFTBDY 

155 

AFTBDY 

156 

AFTBDY 

157 

AFTBDY 

158 

AFTBDY 

159 

AFTBDY 

160 

AFTBDY 

161 

AFTBDY 

162 

AFTBDY 

163 

AFTBDY 

164 

AFTBDY 

165 

AFTBDY 

166 

AFTBDY 

167 

AFTBDY 

168 

AFTBDY 

169 

AFTBDY 

170 

AFTBDY 

171 

AFTBDY 

172 

AFTBDY 

173 

AFTBDY 

174 

AFTBDY 

175 

AFTBDY 

176 



RUN 

CARDS ENCOUNTERED IN INPUr 

UPDATE V?2“^52. 

79/06/22. 

18.29 


KLIMIT - KHXLL 

— ■■ 

AFTBDY 

177 

c+* 



AFTBOY 

178 


INITIALIZE X»Y#2 AND P»Q#R ARRAYS TO 0. 


AFTBDY 

179 




AFTBDY 

160 


DO 10 I ■ l>lLIf1IT 


AFTBDY 

181 


pm • 0*0 


AFTBDY 

182 

10 

CONTINUE 


AFTBDY 

183 


DO 20 J - l^JLIHIT 


AFTBDY 

18<i 


OU) • 0,0 


AFTBDY 

185 

20 

CONTINUE 


AFTBDY 

186 


DO 30 K ■ 1>KLIHIT 


AFTBDY 

187 


R(K> • 0.0 


AFTBOY 

166 

30 

CONTINUE 


AFTBOY 

189 




AFTBDY 

190 


DO 50 I ■ 1>1LIMIT 


AFTBOY 

191 


DO 50 J - I^JLIHIT 


AFTBDY 

192 


DO 50 K • 1, KLIMIT 


AFTBOY 

193 


xmj#K) ■ 0.0 


AFTBDY 

1<)4 


Y(I#J#K) » 0.0 


AFTBDY 

195 


Z(I>J*K) « 0.0 


AFTBDY 

196 

50 

CONTINUE 


AFTBOY 

197 




AFTBDY 

198 


READ PART OF INPUT AND CHECK FOR GRID SIZE. PLOT 

CARDS ARE 

AFTBOY 

199 


READ IN IPLTC ROUTINES. 


AFTBOY 

2 00 




AFTBDY 

201 


READ (5# 1000) IMAX# JMAXf KMAX 


AFTBDY 

202 


WRITE(6>1010) 


AFTBDY 

203 


URITE(6#1020) IMAXi JMAXfKMAX 


AFTBOY 

204 


IF(IMAX.GT.ILIMIT) GO TO 80 


AFTBDY 

205 


IFUMAX.GT.JLIMIT) GO TO 00 


AFTBOY 

206 


IFIKMAX.GT. KLIMIT) GO TO 80 


AFTBDY 

207 


GO TO 100 


AFTBDY 

208 




AFTBDY 

209 

C** 

GRID SIZE specified TOO LARGE 


AFTBOY 

210 




AFTBOY 

211 

60 

CONTINUE 


AFTBOY 

212 


WRITE (6# 1030) IMAX#JHAX#KMAX#ILIMIT# JLIMIT>KLIMIT 


AFTBOY 

213 


STOP 


AFTBOY 

214 

ICO 

CONTINUE 


AFTBDY 

215 




AFTBDY 

216 


REA0(5#I0A0) R1#R2#R3»RA»TERR»ITEPMX 


AFTBDY 

217 


URITE(6>1050) R1«R2>R3*RA«TERR»TTERHX 


AFTBDY 

■ 218 


RE AO (5i 1060) ((X(I»l>l)»Y(I»l>l)*m«l»l))»I>l»INAX) 


AFTBDY 

219 


RE AO (5> 1060) ((X(l«J>l).Y(l»J.l)>Z(l.Jtl))> J-1>JMAX) 


AFTBOY * 

220 




:REATI0M RUN CARDS ENCOUNTERED IN INPUT ■ ■■ ' UPDATE l'.2-452;' 

READ(5#1060) ( ( X( 1^1#K )» Y( 1> lit K) » Z( 1» K ) ) t K-1>KHAX ) ' 

READ(5»1000) IPRT1»IDEL1* JDEH»IPRT2»IDEL2*JDEL2 
WRITE(6/1070) IPRT1»IDEL1»JDEL1»IPRT2*IDEL2»JDEL2 
R£AD(5»10eo» (HED1(I)»I-1,A) 

WRITE(6*1090) IHEDKI). I«l,41 

READ(5>1080) (HED2( I ) • I -1» A) 

WRITEC6>li00) (HED2«n»I-l»A) 

C** 

C*+ CALL INIT TO SET UP INITIAL GUESS AND COMPUTE P»Q,R 

c** 

CALL INIT 

C*» CALL IPRTC TO PRINT INITIAL GUESS IF IPRTl SET. 

C** 

ISOLN - 0 
CALL IPRTC 

C** 

C** CALL IPLTC FOR PLOTTING INITIAL GUESS IF OPTIONS SET. 

CALL IPLTC 

C** COMPUTE COORDINATE SYSTEM 

CALL CALCCR 

C** CALL IPRTC AND IPLTC TO PRINTfPLOT CONVERGED OR 

C^* PARTIALLY CONVERGED SOLUTION IF PROPER OPTIONS SET. 

C4* 

ISOLN ■ 1 
CALL IPRTC 

CALL IPLTC 

C^+ STORE CONVERGED SOLUTION. 

C*+ 

IF(KKK.EQ.O) GO TO ISO 

WRITE(1,1110) ILIMIT# JLIMIT.KLIHIT#IMAXf JMAX#KMAX# 

1 ( (HEDl(n#HED 2 (I))»T«l>A) 

WRITE (1/1120) ((((X(I»J»K)fY(I»J»K)/2(I/JiK))/K«l/KLIHIT)# 

1 J-1#JLIMIT)/I-1/ILIMTT) 

GO TO 200 
C + + 

C++ SOLUTION DID NOT CONVERGE. WRITE A MESSAGE. 

C ++' ‘ . 

150 CONTINUE - 


79/06722. 10. 29. 34. PAGE 6 


AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTfiOY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 
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222 
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22A 

225 

226 
227 
226 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 
240 

249 

250 

251 
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253 

254 

255 

256 

257 
250 

259 

260 
261 
262 

263 

264 
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<1 


CREATION 

RUN 

CARDS ENCOUNTERED IN INPUT 

7 ”’. 79^06/22. 

18.29. 



WRITE(6#1130) YERRMX^IYERIf iyERJ>IYERK»R3#ITERMX 

AFTBOY 

265 




AFTBDY 

266 


C** 

ALL FUNCTION CDMPLeTE 

AFTBDY 

267 



t 

AFTBDY 

268 


200 

CONTINUE 

AFTBDY 

269 



STOP 

AFTBDY 

270 


c#+ 


AFTBDY 

271 



ALL FUNCTIONS COMPLETE. FORMAT STATEMENTS 

AFTBDY 

272 




AFTBDY 

273 


1000 

FGRHAT(16I5) 

AFTBDY 

274 


1010 

F0RMATC«1«) 

AFTBDY 

275 


1020 

F0RMATI»0«>5X#«IMAX#JMAX,KMAX ■ «#3I5) 

AFTBDY 

276 


1030 

F0RMAT(///5X»«*+* GRID SIZE LARGER THAN ALLOWED"# 

AFTBDY 

277 


1 

GRID SIZE INPUT • »»3I5# 

AFTBDY 

278 


2 

/5X>"*** LARGEST SIZE ALLOWED • "»3I5» 

AFTBOY 

279 


3 

/5X»"*f* INCREASE VALUE OF IMXLL* JHXLL.KMXLL IN TEXT"# 

AFTBDY 

280 


4 

« EDITOR”# 

AFTBDY 

281 


5 

/5X»"*** RUN ABORTED ") 

AFTBDY 

282 


10^0 

F0RMAT(AF10.0»2I5) 

AFTBOY 

283 


1050 

F0RMAT(”0”#5X#”Rl#R2#R3#R4f I£RR#ITFRMX ■ ”#4E15.5» 215) 

AFTBOY 

284 


1060 

FQRHATOFIO.O) 

AFTBDY 

285 

N) 

|_j 

1070 

F0RMAT(«0»#5X#”IPRTl»I0ELlf JDEL1»IPRT?#IDEL2# JDEL2 ■ ”#615) 

AFTBDY 

286 

1060 

F0RHAT(4A10) 

AFTBDY 

287 


1090 

FORMAT(»0”i5X#”HEDl ■ ”f4A10) 

AFTBDY 

288 


1100 

F0RMAT(»0«#5X#”HED2 - »#4A10) 

AFTBDY 

289 


1110 

F0RMATC6I5#8A10) 

AFTBOY 

290 


1120 

FQRMAT(1X,6E15.6) 

AFTBOY 

291 


1130 

F0RMAT(»1”////5X#»**^ SOLUTION DID NOT CONVERGE. SOLUTION”# 

AFTBDY 

292 


1 

» NOT STORED ON UNIT 1”//5X#”*^^ MAX ERROR IN Y ■ ”#E15.6# 

AFTBOY 

293 


2 

••DCCUREO AT LOCATION « »# I5#«#«» r5»»#«#I5//5Xf «♦♦♦ MAX«# 

AFTBOY 

294 


3 

" ERROR ALLOWED ■ "lE 15.6//5X#«*<'^ NO. OF ITERATIONS”# 

AFTBDY 

295 



” PERFORMED - ”#I5) 

AFTBOY 

296 


c** 


AFTBDY 

297 



END 

AFTBOY 

298 



SUBROUTINE INIT 

AFTBDY 

299 


C** 


AFTBDY 

300 



THIS SUBROUTINE COMPUTES INITIAL GUESS ON THE AFTER BODY 

AFTBDY 

301 



MODEL. ONLY ONE-HALF OF THE BODY IS CONSIDERED. THE BODY 

AFTBOY 

302 



CONSISTS OF A 4-DEGREE WEDGE SECTION AND A AFTER BODY 

AFTBOY 

303 



WHICH IS CURVED IN X AS WELL AS Z-DIRFCTION. THE WEDGE 

AFTBDY 

304 



HAS A LENGTH OF 62.3 AND THE AFT BODY HAS A LENGTH OF ’ 

AFTBOY 

’305 



27.7. THE HALF WIDTH 13 12. IN THE AFTER BODY 

AFTBOY 

306 



Y-CENTERLINE (YCL) IS GIVEN BY 

AFTBDY 

307 



YCL ■ C > ID^O - CX-E)^CX-E))**0.5 

AFTBDY 

308 
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REATION RUN 

C** 

C** 

C** 

C** 

C** 

C+* 

€♦♦ 

C** 

C** 

C++ 


c*t 

c*+ 


20 


AO 

50 


C*+ 
C++ 
C + + 
C++ 


75 


CAROS ENCOUNTERED IN INPUT ' 

WHERE C - -A57.8350 

0 - A63.316A " 

E ■ 0A.5670A 

IN LATERAL DIRECTION 

V ■ YCL - R + (R+R - Z+7) +♦ 0,5 
WHERE R ■ (A+B(Y-62.3)) ++-1.0 

A - 3.5521E-A P-2.23505E-A 

NOTEi Z COMPUTED FROM CENTERLINE OF BODY. 

COMMON /0XXX/ XIIMXLL* JMXLL»KMXLL),Y(IHXLL, JMXLL.KMXLU* 

1 ZdMXLL* JMXLL,KMXLL).P(IMXLL)»0( JMXLU»R(KMXLL) 

COMMON /AXXX/ IMAX, JMAX , KMAX , IL IM I T» JL IMIT»KL I MI T# ISOLN# 

1 KKK,XERRMX,YERRMX»ZERRMX,IXERI.IXERJdXERK,IYERI,IYERJ» 

2 lYERK# IZERI,IZERJ» TZERK>R1,R2»P3»RA, IERR» ITERMX, IPRTl# 

3 I0EL1» JDELl# IPRT2» I0EL2. J0EL2»NPLT1» ITYPE, ISTN,SIZE1» 

A SIZE2,LL>HE01l A1,HE02( A) ,IPLSET,FILLXXJ25» 

COMPUTE YCL. TANA IS TANGENT PF A-DEGREFS. 

TANA - .0699268119 
A ■ 0.00035521 
B • 0.000223505 
DO 50 I - l.IMAX 
Z(I.l.l) - 0.0 

IF(X(I.1>1).GT.62.3) GO TO 20 
Y(l,ld) • X(I»1.U ♦ TANA 
GO TO AO 
CONTINUE 

Y(I»1/1» • -A57.8350 ♦ (A63.316A+A63.316A - 
L (X(l4l>ll - 9A.5670AI+ (X(Itl.l) - 9A.5670AM ++0.5 

CONTINUE 
CONTINUE 

ON ALL I STATIONS X IS CONSTANT AND ON ALL 
K-STATIONS Z IS CONSTANT 

DO 75 I • l.IMAX 
DO 75 J ■ 1»JMAX 
DO 75 K ■ l.KMAX 
XII, J.K) - Xlld.l) 

Z(I>J,K) ■ Z(l.l.K) 

CONTINUE 


UPDATE 1.2-A52V 


c** 


79/06/22/ 

16.29.3A 

AFTBDY 

309 

AFTBDY 

310 

AFTBDY 

311 

AFTBDY 

312 

AFTBDY 

313 

AFTBDY 

31A 

AFTBDY 

315 

AFTBDY 

316 

AFTBDY 

317 

AFTBDY 

310 

AFTBDY 

319 

AFTBDY 

320 

AFTBDY 

321 

AFTBDY 

322 

AFTBDY 

323 

AFTBDY 

324 

AFTBDY 

325 

AFTBDY 

326 

AFTBDY 

327 

AFTBDY 

328 

AFTBDY 

329 

AFTBDY 

330 

AFTBDY 

331 

AFTBDY 

332 

AFTBDY 

333 

AFTBDY 

334 

AFTBDY 

335 

AFTBDY 

336 

AFTBDY 

337 

AFTBDY 

338 

AFTBDY 

339 

AFTBDY 

340 

AFTBDY 

341 

AFTBDY 

342 

AFTBDY 

343 

AFTBDY 

344 

AFTBDY 

345 

AFTBDY 

346 

AFTBDY 

347 

AFTBDY 

346 

AFTBDY 

349 

AFTBDY 

350 

AFTBDY 

351 

AFTBDY 

352 
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c c 


CREATION RUN 


CAROS ENCOUNTEREn IN INPUT " '“ UPDATE 1. 


C*^ ON J • JMAX V • YC1»JMAX»1I ^ ^ 

C** 

DO TOO I - l^IHAX 
DO 100 K • 1#KMAX 

Y(I,JMAX#K) ■ Y(1#JMAX»1) 

100 CONTINUE 

C^^ COMPUTE Y ON J • 1. UPTO X • 62.3 Y STAYS CONSTANT IN K 

C+^ DIRECTION. FOR X GREATER THAN 62.3 Y HAS TO BE COMPUTED. 

C** ^ 

DO 130 I • I^IMAX 
DO 120 K - 2»KMAX 

IF{X(I#1>1).GT.62.3) 60 TO 110 
Y(I#1#K) ■ YII#1#1) 

GO TO 115 
110 CONTINUE 

RR • 1.0 / (A ♦ B^(X(Ifl»n - 62.3)1 
Y(I#1#K) • Yn»l#l) - RR ♦ (RR + RR - 2(I»1#K) 

1 ♦ 2(I#1#K)) ♦♦ 0.5 

115 CONTINUE 
120 CONTINUE 
130 CONTINUE 
€♦♦ 

C^* COMPUTE Y FROM J • 2 THRU JMAX - 1. POINTS ARE 

C*^ DISTRIBUTED IN THE SAME RATIO AS ON I • 1#K-1 LINE. 

C^* 

J1 • JMAX - 1 
DO 150 I - I^IMAX 
DO 150 J • 2»J1 
DO *150 K « 1»KMAX 

YCI#J»K) - Y(I#1»K) ^ (YlIfJHAXfK) - Y(I»1#K)) ♦ 

1 (Y(1>J#1) * Y(l»l»in /(Y(1#JMAX#1) - Y(l»l#l)) 

150 CONTINUE 

€♦♦ INITIAL GUESS SPECIFIED. COMPIJTF P,Q, AND R 

C** 

II ■ IMAX - 1 
DO 200 I - 2#I1 

xuu - xiuifipi) - 2.04<x(i#i»i) ♦ xn-ifiti) 

XU • 0.5 ♦ CX(I+1>1»1) - X(I-l#l#D) 

P(I) • - xuu / (XU^XU^XU) 

200 CONTINUE 

DO 250 J • 2# J1 


AFTBDY 353 
AFTBOY 35A 
AFTBDY 355 
AFTBOY 356 
AFTBDY 357 
AFTBOY 350 
AFTBOY 359 
AFTBOY 360 
AFTBDY 361 
AFTBOY 362 
AFTBDY 363 
AFTBDY 364 
AFTBOY 365 
AFTBOY 366 
AFTBDY 367 
AFTBDY 368 
AFTBOY 369 
AFTBOY 370 

AFTBDY 371 

AFTBDY 372 
AFTBDY 373 
AFTBOY 374 
AFTBDY 375 
AFTBOY 376 
AFTBDY 377 
AFTBDY 378 
AFTBOY 379 
AFTBOY 380 
AFTBOY 381 
AFTBDY 302 
AFTBOY 303 
AFTBDY 384 
AFTBDY J 385 
AFTBDY 386 
AFTBDY 387 
AFTBOY 388 
AFTBOY 389 
AFTBDY 390 
AFTBDY 391 
AFTBDY 392 
AFTBDY 393 
AFTBOY 394 
AFTBDY ■ 395 
AFTBOY 396 


ATION RUN CARDS ENCOUNTERED IN INPUT UPDATE 1.2-452 

yUU ■ Y(1»J+1.1) - 2.0*Y(l» J.l) ♦ Y(1»J-1,1) 

YU • 0.5 * (Y(lfJ4^1>ll - YCl«J-l»in 

0(J) » - YUU / (YU*YU*YU> 

250 CONTINUE 
C** 

KI ■ KMAX - 1 • 

DO aOO K ■ 2>K1 

ZUU » 2(1,1,K+1» - 2.0*m,l,K> + 2(l»l,K-l) 

2U • 0.5 * (2(1»1.K+1» - Z(1,1,K-1)) 

R<K) • - 2UU / (2U*2U*ZUI 
300 CONTINUE 
C** 

RETURN 
END 

SUBROUTINE CALCOR 
C** 

C** THIS SUBROUTINE COMPUTES Y-CORDINATE ONLY IN THE INNER GRID. 

C** X AND 2 COORDINATES IN THE GRID ARE FIXED. FOR SIMPLICITY 

C** U SIGNIFIES XI» V - ETA» AND W - TAU. THUS XU SIGNIFIES 

C** DX/OXI. XUU • D2X/DXI2 AND SO ON. 

C** KKK ■ 0 SIGNIFIES UNCONVERGED SOLUTION. 

C** 

COMMON /BXXX/ X ( IMXLL* JMXLL»KMXLL ) » Y( IMXLL> JMXLL,KMXLL ) » 

1 2(IMXLLi JMXLL«KMXLL)»P(IMXLH>0«JMXLL»»R(KHXLL) 

COMMON /AXXX/ IMAX, JMAX, KMAX , IL I MI T, JL IMI T» KL IHI T» ISQLN, 

1 KKK.XERRMX, YERRMX.ZERRMX, IXERI.TXERJ»IXERK>IYERI>IYERJ» 

2 IYERK>I2ERI»I2ERJ»IZERK»R1»R2.R3.RA»IERR» ITERMX* IPRTl# 

3 lOELl. JDEH.IPRT2#IDEL2» J0EL2.NPLT1,ITYPE»ISTN.SI2E1* 

4 SI2£2»LL/HED1(4»,HED2(4>, IPLSET,FILLXX(25) 

C*+ 

ISOLM ■ 1 

KKK - 0 
lERFLG - 0 
II ■ IMAX - 1 

J1 • JMAX - 1 

K1 - KMAX - 1 

DO 600 UK • 1»ITERMX 
LL - UK 
YERRMX ■ 0.0 

DO 400 I ■ 2«I1 

DO 400 J ■ 2,J1 

DO 400 K ■ 2»K1 

YTEMP - Y(I.J#K) 

PPP • P(I) 


t 


79/06/22* 18,29,34. PAGE 10 

AFTBDY 397 

AFTBOY 398 

AFTBDY 399 

AFTBDY 400 

AFTBDY 401 

AFTBOY 402 

AFTBDY 403 

AFTBDY 404 

AFTBDY 405 

AFTBDY 406 

AFTBOY 407 

AFTBOY 408 

AFTBDY 409 

AFTBOY 410 

AFTBDY 411 

AFTBDY 412 

AFTBOY 413 

AFTBDY 414 

AFTBDY 415 

AFTBOY 416 

AFTBDY 417 

AFTBDY 418 

AFTBOY 419 

AFTBDY 420 

AFTBDY 421 

AFTBOY 422 

AFTBDY 423 

AFTBDY 424 

AFTBOY 425 

AFTBDY 426 

AFTBDY 427 

AFTBDY 428 

AFTBDY 429 

AFTBDY 430 

AFTBOY 431 

AFTBDY 432 * 

AFTBOY 433 

AFTBDY 434 

AFTBDY 435 

AFTBDY 436 

AFTBDY 437 

AFTBOY 438 

AFTBDY 439 

AFTBDY 440 


t) 


I 







* in 


CREATION RUN 


CARDS ENCOUNTERED IN INPUT ^ ^ ^ 18.29.3A. * / PAG 


1 

1 

1 

C^* 

1 

1 

1 

1 

1 

1 

1 

1 


QQQ 

■ 0( J) 

* ’ ' " • - ■ - 

AFTBDY 

4A1 

RRR 

• R(K) 


AFTBOY 

442 

IPl 

- I + l 


AFTBDY 

443 

IMl 

- I-l 


AFTBDY 

444 

JPl 

• J+1 


AFTBDY 

445 

JMl 

• J-1 


AFTBDY 

446 

KPl 

■ K*l 


AFTBDY 

447 

KMl 

■ K-1 


AFTBDY 

448 

XU 

(X(1P1»J,K) - X(IMl,J,Kn ♦ 0,5 


AFTBDY 

449 

YU 

(Y(IPl»J>K) - Y(IH1>J>K)) * 0.5 


AFTBDY 

- 450 

ZU 

(ZdPlfJfK) - Z(TMlfJ,Kn ♦ 0.5 


AFTBDY 

451 

XV 

(X(I>JP1,K) - XII*JH1#K)) ♦ 0.5 


AFTBDY 

452 

YV 

CY<dJPl,K) - Y(I*JHl»Kn ♦ 0.5 


AFTBDY 

453 

2V 

(Zd^JPl^K) - Z(IfJMl^K)) ♦ 0.5 ■ 


AFTBDY 

454 

XU 

(Xd,J,KPl) - Xd#J#KHD) ♦ 0.5 ^ 


AFTBDY 

455 

YW 

(Y(I»J»KP1) - Y(I,J,KMD) ♦ 0.5 


AFTBDY 

456 

zw 

(Z(I>J»KP1) - ZCr.J.KMD) * 0.5 


AFTBDY 

457 

YUU 

• Y«IP1»J»K) - 2.0*Y(I,J,K) ♦ Y(IM1,J,K) 


AFTBDY 

458 

YVV 

■ Yd, JPl, K) - 2.0^Y(I,J#K) ♦ YCI,JM1,K) 


AFTBDY 

459 

yuw 

• Y(I>J»KP1» - 2.0*Y(I,J,K) + YJI»J»KH1) 


AFTBDY 

460 

YUV 

• 0.25*(Y(IP1. JPl.K) - YtTPl.JHl^KJ - YIIM1» 

JP1»KI + 

AFTBDY 

461 


Y(IMl.JHl.K)) 


AFTBOY 

462 

YUW 

- 0.25*(Y(IP1» J»KP1) - Y(IP1»J»KH1I - Y(IMI»J»KP1) ♦ 

AFTBOY 

463 


YdHl,J,KMD) 


AFTBDY 

464 

YVW 

« 0.25*(Yd,JPl,KPl) - Y(I,JH1,KP1) - Y(I,JP1,KH1) + 

AFTBOY 

465 


Yd,JMl,KMD) 


AFTBOY 

466 




AFTBDY 

467 

AJ 

■ XU*(YV*ZW - ZV^YW) - XV*(YU*ZW - YW*2U) 4- 


AFTBDY 

468 


XW*(YU+ZV - YV*ZU) 


AFTBOY 

469 

CUU 

• <XV*YW - XW*YV)**2 ♦ CXVAZW - XW*ZV)^*2 + 


AFTBDY 

470 


(YV*ZW - YW*ZV)**2 


AFTBOY 

471 

CUV 

■ (XU*YW - XW*YU»*(XV*YW - XW+YV) + (XU*ZW - 

XW*ZU) ♦ 

AFTBDY 

472 


(XV*ZW - XW+ZV) + (YU*ZW - YW*ZU) ♦ (YV+ZW - 

YW+ZV) 

AFTBOY 

473 

cuw 

• (XU*YV - XV*YU)*{XV4‘YW - XW*YV) * CXU*2V - 

XV^ZU) ♦ 

AFTBDY 

474 


CXV*ZW - XW^ZV) f (YU+ZV - YV+ZU) ♦ lYV+ZW - 

YW*ZV) 

AFTBDY 

475 

cvv 

- (XU*YW - XW*YU)**2 ♦ (XU+7W - XW*ZU»+*2 ♦ 


AFTBDY 

476 


(YU+ZW - YW*7U)**2 


AFTBDY 

477 

cvu 

« (XU*YV - XV+YU)*(XU*YM - XW*YU»* (XU+ZV - 

XV+ZU) + 

AFTBDY 

478 


(XU*ZW - XW*ZU) ♦ (YU+ZV - YV*ZU> * (YU*ZW 

- YW+ZU) 

AFTBOY 

479 

CUU 

• (XU^YV - XV*YU)**2 4- (XU^ZV - XV*ZU)**2 + 


AFTBOY 

4 80 


<YU+ZV - YV+ZU)++2 


AFTBOY 

481 

CF 

• 0.5 / (CUU ^ CVV + CWW) 


AFTBDY 

482 

YTEMPl « CF+<CUU+«Y(Uli J,K) + Y«I-1,J,K)) - 2.0*CUV+YUV 

AFTBOY 

483 


♦ 2.0 + CUW+YUW ♦ CVV*(Yd, J + 1>K) + Yd,J-l,K)) 


AFTBDY 

484 







11 T” 



:ATiaN RUN 


CAROS ENCnUNTERED IN INPUT UPDATE l.2-^52. 79/06/22. 18. 29.3A. 


2 ♦ CWU*(Y(I*J»K+1) ♦ Y(I«J>K-in - 2.0fCVW*YVW ♦ 

3 (PPP*YU + QQQ+YV + RRR+YWJ ♦ (AJ*AJ)I 
Y(I,J,K) • (1.0-Rl) ♦ YTEHP + R1*YTEMP1 

COMPUTE MAX Y-ERROR 

YERR - ABS(YTEMP - Y(I.J.K)) 

IF(YERR.LE.YERRMX) GO TO 150 
YERRMX ■ YERR 
lYERI • I 
lYERJ • J 
lYERK - K 
CONTINUE 
CONTINUE 

CHECK IF ERROR NEEDS TO BE PRINTED 

lERFLG - lERFLG ♦ 1 
IF(IERFLG.LT.IERR) GO TO A50 

WRITE (6.1000) LL. YERRMX, IYER I. IYER J.IYERK 
lERFlG - 0 
CONTINUE 

CHECK FOR CONVERGENCE 

IF( YERRMX. LT.R3) KKK - 1 
IF(YERRMX.LT.R3) GO TO 650 
600 CONTINUE 
C + * 

650 CONTINU'E 

C** 

RETURN 

1000 F0RMAT("0",5X, "ITERATION NO., MAX Y-ERR, LOC • ", 15, 2X, E15.5, 2X, 
1 316) 

C+ + 

END 

SUBROUTINE IPRTC 
C** 

C** SUBROUTINE USED TO PRINT INITIAL, CONVERGED OR PARTIALLY 

C** CONVERGED SOLUTION. ALL PERTINENT PARAMETERS ARE PASSED 

C*+ THRU COMMON. 

ISQLN - D INDICATES INITIAL GUESS 

C**' ISOLN » 1 AND KKK . 0 INDICATES PARTIALLY CONVERGED SOLN 

C** ISOLN • 1 AND KKK - 1 INDICATES CONVERGED SOLUTION. 


C** 

C** 

C** 


150 
A 00 
C** 


450 

C++ 

C* 


AFTBOY 

AFTBDY 

AFTBOY 

AFTBOY 

AFTBOY 

AFTBDY 

AFTBOY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBOY 

AFTBDY 

AFTBOY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBOY 

AFTBOY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBOY 

AFTBOY 

AFTBOY 

AFTBDY 

AFTBOY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBDY 

AFTBOY 

AFTBDY 

AFTBOY 

AFTBDY 

AFTBOY 

AFTBOY 

AFTBDY 

AFTBOY 


485 

466 

487 

488 
4B<) 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

511 

512 

513 

514 

515 

516 

517 

518 

519 

520 

521 

522 

523 

524 

525 

526 

527 
526 
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I) 


I 
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CREATION 

RUN 

CARDS eNCnUNTERED IN INPUT / “ ""Zll ... 

79/06/22. 

16 •2.*’ *34 


c** 


AFTBDY 

529 

... 


COHHON /BXXX/ X ( IHXLL> JHXLL«KMXU ) * Y( IMXLL* JHXLL«KHXLU * 

AFTBOY 

530 


1 

Z(IHXLL> JHXLL»KHXLL)*PIIHXLU>0( JHXLL)fR(KHXLL) 

AFTBOY 

531 



COHHON /AXXX/ IHAX« JHAX.KHAX« ILIHITt JLIHIT»KLIHIT> ISOLN» 

AFTBOY 

532 


1 

KKK>XERRHX>YERRHX>7ERRHX>IXERI*TXFRJ*IXERK>IYERI*IYERJ* 

AFTBDY 

533 


2 

IYERK>IZERI#IZERJ»rZERK#Rl#R2»R3>PA,lERR#ITERMX#IPRTl» 

AFTBOY 

534 


3 

IDELl^JD£Ll#IPRT2^IDEL2,J0EL2#NPLTl,tTYPE#ISTN#SIZEl> 

AFTBDY 

535 


A 

SI 2E2#LL#HED1C^)#HE0R< A I PLSE T, F ILL XX (2 5 ) 

AFTBDY 

536 




AFTBDY 

537 



IF(ISOLN.EQ.l) GO TO 10 

AFTBDY 

538 



IPRTX - IPRTl 

AFTBDY 

539 



lOELX « IDELl 

AFTBDY 

540 



JDELX - JDEH 

AFTBOY 

541 . 



GO TO 20 

AFTBDY 

542 


10 

CONTINUE 

AFTBDY 

543 




AFTBOY 

544 



IPRTX • IPRT2 

AFTBOY 

545 



IDELX - IDEL2 

AFTBDY 

546 



JDELX • JDEL2 

AFTBOY 

547 


20 

CONTINUE 

AFTBDY 

548 




AFTBDY 

549 

o 


IFUPRTX.EQ.O) GO TO 600 

AFTBDY 

550 




AFTBDY 

551 



WRITE HEADINGS 

AFTBDY 

5 52 




AFTBDY 

553 



URITE(6#1000) HEOl 

AFTBDY 

554 



WRITE(6f ICIO) HE02 

AFTBDY 

555 



IF(ISOLN.EO.O) WRITE<6#1020) 

AFTBOY 

556 



IF(ISOLN.EQ.l.ANO.KKK.EO.O) WR I TE ( 1030 ) 

AFTBOY 

5 57 



IFdSOLN.EQ.l.AND.KKK.EO.n WRITE ( 6# 10A01 

AFTBDY 

558 

..... 


WRITEC6#1050) IMAX, JMAX# KMAX 

AFTBOY 

559 



IF(ISOLN.EQ.O) GO TO 50 

AFTBOY 

560 



WRITE(6#1060) LL 

AFTBDY 

561 



WRITEC6d070) YERRMX, IYER I » IYER J, lYERK 

AFTBDY 

562 


50 

CONTINUE 

AFTBDY 

563 


c*+ 


AFTBDY 

■ 564 



DO 300 I « IMHAX/ IDELX 

AFTBOY 

565 



DO 290 J • 1#JHAX>J0ELX 

AFTBDY 

566 



WRITE (6>10B0) IfJ 

AFTBOY 

567 



WRITE(6«1090) ((X(I«J>KI»Y(T»J»K),Z(I»J,Kn>K-l*KHAX) 

AFTBDY 

568 


290 

CONTINUE ' 

AFTBOY 


. ... 

300 

CONTINUE 

AFTBOY 

570 




AFTBOY 

571 



600 

CONTINUE " ' * 

AFTBOY "■ 

572 “ 




AT ION RUN CARDS ENCOUNTERED IN jnPUT UPDATE ■1>2-A52V“""~ 79/06/22^Ib;29.3 


RETURN 

C** 

C** FORHAT STATEMENTS ' ' 

C*+ 

1000 FCRMATC'1"//A0X»AA10) 

1010 F0RMAT(/A0X>AA10) 

1020 F0RMATC'0''>A9X»''*** INITIAL SOLUTION *♦*") 

1030 F0RMAT("0'',A1X#«*** PARTIALLY CONVERGED SOLUTION ♦♦♦'•) 
lOAO F0RMAT("0">A8X>''*** CONVERGED SOLUTION ♦♦♦••» 

1050 F0RMATC'0'',A1X>''IMAX • '•> 13, AX, •' JMAX ■ 13, AX,*'KHAX ■ «,13) 

1060 FORMATC'O", ATX, "ITERATION PERFORMED • ",15J 

1070 FDRHATC ",3AX,"MAX Y-ERR ■ ",EIA,6," AT LOC I,J,K ■ ",3I5) 
1U80 FORMAT t //5X,"1 • »,I3," J ■ ",I3," (X,Y,I)") 

1090 F0RMAT(5X,E11.5,3X,E11,5,3X,E11.5,3X,E11.5,3X,E11.5,3X,E11.5, 
1 3X,E11.5,3X,E11.5,3X,E11.5) 

C** 

END 

SUBROUTINE IPLTC 

C** 

C++ SUBROUTINE USED TO PLOT INITIAL CONVERGED OR PARTIALL 

C*» CONVERGED SOLUTION, SEE COMMENTS UNDER IPRTC ROUTINE 

c *♦ 

COMMON /BXXX/ X ( IMXLL, JMXLL ,KMXLL » , YHMXLL, JMXLL,KMXLL ) , 

1 Z(IMXLL,JMXLL,KMXLU,P(IMXLL»,0(JMXLL),R(KMXLU 
COMMON /AXXX/ IMAX, JMAX , KM AX, I L IM IT, JL IM IT,KL IMIT, ISOLN, 

1 kkk,xerrmx,yerrmx,zerrmx,ixeri,ixerj,ixerk,iyeri,iyerj, 

2 lYERK, IZERl, IZERJ, I ZERK,R1 ,R2, R 3,R A, TERR, I TERMX, IPRTl, 

3 lOELljJDELl, IPRT2, IDEL2, JDEL2,NPLT1, ITYPE, ISTN,SIZE1, 

A SIZE2,LL,HEDl(A),HEn2(A),IPLSET,FILLXXI25> 

C** 

DIMENSION X0UM<75), YDUM(75),Z0UM(7F),HE0( A) 

1 

C** 

C+* READ A CARD AND SEE IF PLOTS PLOTTED OR NOT 

C*+ 

READ(5,1000) NPLTl 
IF (NPLTI.EO.O) go to 900 

C** 

IF(IPLSET.EQ.l) 60 TO 20 
IPLSET « 1 
CALL PSEUDO 

CALL FONTS (1) - 

20 CONTINUE 

C** 


AFTBOY 573 
AFTBDY 57A 
AFTBDY 575 
AFTBOY 576 
AFTBOY 577 
AFTBDY 578 
AFTBDY 579 
AFTBDY 580 
AFTBOY 581 
AFTBDY 582 
AFTBDY 583 
AFTBDY 58A 
AFTBDY 585 
AFTBDY 586 
AFTBDY 587 
AFTBDY 588 
AFTBOY 589 
AFTBOY 590 
AFTBDY 591 
AFTBDY 592 
AFTBDY 593 
AFTBOY 59A 
AFTBDY 595 
AFTBDY 596 
AFTBOY 597 
AFTBDY 598 
AFTBDY 599 
AFTBDY 600 
AFTBDY 601 
AFTBOY 602 
AFTBOY 603 
AFTBDY 60A 
AFTBOY 605 
AFTBDY 606 
AFTBOY 607 
AFTBDY 608 
AFTBDY 609 
AFTBDY 610 
AFTBOY 611 
AFTBDY 612 
AFTBDY 613 
AFTBOY 614 
AFTBOY 615 
AFTBDY 616 




J) 


I 





CREATION RUN 
C** 




CARDS ENCOUNTERED IN INPUT ' 


UPDATE !• 2-^52.’ 


‘79/06/22.^18. 29.3A 




READ NPLTl CARDS AND PLOT 

DO 500 KLH • 1/NPLTl 

READC5#10fi0) ITYPE#rSTN»SIZEl 

PRINT HEADINGS. 

CALL SCRIBE(0.3f 1.5>0.15»0.5>HE01#90.»40#9) 
CALL SCRIBE(0.55#1.5>0.l5#0.5*HED2t90.^A0»9) 

IFCISOLN.EQ.l) GO TO 50 

H£D(1) ■ lOH INITIA 

HED(2) ■ lOHL GUESS " 

HEDJ3) • IH- 
INUHC -21 
GO TO 100 




5C 


K5 

VO 


CONTINUE 

IF(KKK.EQ.l) GO TO 75 
HEO(U 
HEDC2) 

HED(3) 

HEO<^) 

INUMC 


lUH PARTIA 

lOHLLY CONVER 
lOHGED SOLUTI 

lOHON 

37 


€♦>«> 


75 


C*t 

100 

C** 

C# + 
C** 


GO TO lOO 
CONTINUE 

HEOCl) ■ lOH CONVER 

HEDC2) ■ lOHGED SOLUTI 

HE0I3) • 7H0N 

INUMC • 27 

CONTINUE 

CALL SCRIBE(0.70f3.0#0.1#0.5»HEO#9D.O>INUHC#9) 
WRITE ARRAY SIZES 

ENCODE (A0^1020>HED) IMAX» JMAX»KHAX 
CALL SCRIBE(0.85#3.0#0.1#0.5fHED»90.0f36»9l 
IFnTYPE.EQ.2) GO TO 1^0 
XX • X(ISTN#1#1) 

ENCGDE(A0»1060>HE0) ISTN>XX 
GO TO 160 




AFTBDY 

617 

AFTBDY 

618 

AFTBDY 

619 

AFTBDY 

620 

AFTBDY 

621 

AFTBDY 

622 

AFTBDY 

623 

AFTBDY 

624 

AFTBDY 

625 

AFTBDY 

626 

AFTBDY 

627 

AFTBDY 

626 

AFTBDY 

629 

AFTBDY 

630 

AFTBDY 

631 

AFTBDY 

632 

AFTBDY 

633 

AFTBDY 

634 

AFTBDY 

635 

AFTBDY 

636 

AFTBDY 

637 

AFTBDY 

636 

AFTBDY 

639 

AFTBDY 

640 

AFTBDY 

641 

AFTBDY 

642 

AFTBDY 

643 

AFTBDY 

644 

AFTBDY 

645 

AFTBDY 

646 

AFTBDY 

647 

AFTBDY 

646 

AFTBDY 

649 

AFTBDY 

650 

AFTBDY 

651 

AFTBDY 

652 

AFTBDY 

653 

AFTBDY 

654 

AFTBDY 

655 

AFTBDY 

656 

AFTBDY 

657 

AFTBDY 

656 

AFTBDY 

659 

AFTBDY “ 

660 



KEATION RUN 

no 


IbC 


C** 

C** 

€♦♦ 

C** 

C** 

C** 

C** 


180 

C** 

C** 

C+* 


190 


1 

200 


C+* 

C*+ 

C++ 


CARDS ENCOUNTERED IN INPUT ; UPDATE 1.2-A52'. 


CONTINUE - 

n ■ Z(1/1»ISTNI 

ENC0DEU0»10e0»HED) ISTN»ZZ 

CONTINUE 

CALL SCRIBEa.2«3.0>0,l«0.5*HED>90.0«38*9) 

CALL NFRAME 

START PLOT. IF ITYPE • 1 THEN Y-Z PLANE PLOTTED. IF ITYPE -2 
THEN X-Y PLANE PLOTTED. 

IF(ITYPE.EQ.2) GO TO 3D0 

COMPUTE YMIN*YMAX» ZMIN»ZMAX 

YMIN ■ 10000.0 
YMAX « -10000.0 
ZHIN ■ 10000.0 
ZHAX • -10000.0 
DO 160 J ■ 1>JMAX 
DO 180 K ■ IfKHAX 

IF(Y(ISTN.J>K).LT.YMIN) YMIN • Y(ISTN>J»K) 

1F( YdSTN. JiKI.GT.YMAX) YMAX • Y(tSTN»J»K» 

IF(Z(1STN» J»K).GT.ZMAXI ZMAX • 7«ISTN»J#KI 
IF<Z(1STN/J»K).LT.ZMIN» ZMIN » ZMSTN»J»K» 

CONTINUE 

SIZE2 ■ SIZEl + (YMAX - YMIN) / (ZMAX - ZHIN) 

START PLOT. ETA LINES FIRST, 

DO 200 J « 1,JMAX 
DO 190 K • 1»KMAX 

YDUM(K) • Y(ISTN,J»K) 

ZDUM(K) ■ ZdSTN.J.K) 

CONTINUE 

CALL INF0PLT(0,KMAX,ZDUM,1.Y0UM,1»ZMIN»ZMAX,YMIN»YMAX, 
1.0>0»HE0»0>HED,0»SIZEl>SIZE2»0.0,2.O) 

CONTINUE 

ZETA LINES. 


KK • 0 

DO 230 K ■ 1»KMAX 

IF (K.EQ.KMAX) KK • 1 
DO 216 J ■ l.JMAX 


79/06/22, 

ie«29«3^ 

AFTBOY 

661 

AFTBDY 

662 

AFTBDY 

663 

AFTBDY 

664 

AFTBOY 

665 

AFTBDY 

666 

AFTBDY 

667 

AFTBDY 

668 

AFTBDY 

669 

AFTBOY 

670 

AFTBDY 

671 

AFTBDY 

672 

AFTBDY 

673 

AFTBOY 

674 

AFTBDY 

675 

AFTBDY 

676 

AFTBDY 

677 

AFTBOY 

678 

AFTBDY 

679 

AFTBOY 

680 

AFTBDY 

681 

AFTBDY 

682 

AFTBOY 

683 

AFTBOY 

684 ' 

AFTBOY 

685 

AFTBOY 

686 

AFTBDY 

687 

AFTBOY 

688 

AFTBOY 

689 

AFTBOY 

690 

AFTBOY 

691 

AFTBOY 

692 

AFTBOY 

693 

AFTBDY 

694 

AFTBDY 

695 

AFTBOY 

696 

AFTBDY 

697 

AFTBOY 

698 

AFTBOY 

699 

AFTBOY 

700 

AFTBOY 

701 

AFTBDY 

702 

AFTBDY 

703 

AFTBOY 

704 


«l I 


li 


I 



CREATION RUN 

CAROS ENCOUNTERED IN INPUT ' J “UPDATE ' 

"79/06/22. 

18.29- 



YDUMCJ) • YCISTN» J#K) 

AFTBDY 

705 



ZOUN(J) - Z(ISTN>J>KI 

AFTBDY 

706 


215 

CONTINUE 

AFTBOY 

707 



CALL INFCPLT(KK»JHAX>ZDUN*l»yOUH>l>ZHIN»ZHAX«yMIN*YHAX» 

AFTBDY 

708 


1 

1.0«0«HEO*0/HED>0«SIZE1>SIZE2»0.0»2.0) 

AFTBDY 

709 


230 

CONTINUE 

AFTBDY 

710 



GO TO 480 

AFTBDY 

711 




AFTBDY 

712 



' • 

AFTBDY 

713 



X-Y PLANE IS TO BE PLOTTED. 

AFTBDY 

714 




AFTBOY 

715 


300 

CONTINUE 

AFTBOY 

716 



XHIN ■ 10000.0 

AFTBOY 

717 



XHAX - -10000.0 

AFTBDY 

718 




YMIN • 10000.0 

AFTBDY 

719 



YMAX • -lOCOO.O 

AFTBOY 

720 



DO 320 I - l^IMAX 

AFTBOY 

721 



DO 320 J ■ 1»JMAX 

AFTBDY 

722 

.. 


IF(xn>J>ISTN».LT.XHIN) XMIN • X(I»J>ISTN» 

AFTBOY 

723 



IF(X(I/J>ISTN).GT.XMAX) XHAX • X(T#J»tSTN) 

AFTBOY 

724 



IF(ya>J>ISTN».LT.yMIN) YMIN • y(I.J»ISTN» 

AFTBOY 

725 



IF(Y(I> JfISTNl.6T.yMAX) YMAX - Yd.JflSTNI 

AFTBDY 

726 


320 

CONTINUE 

AFTBOY 

727 



SIZE2 ■ SIZEl * (YMAX - YHIN) / IXHAX - XHIN) 

AFTBOY 

728 




AFTBDY 

729 



ETA LINES 

AFTBOY 

730 




AFTBDY 

731 



DO 350 J • IfJMAX 

AFTBOY 

732 



DO 340 I ■ IfIMAX 

AFTBDY 

733 



XOUM(I) - XdfJfISTN) 

AFTBDY 

734 



YOUMd) ■ YdfJ.ISTN) 

AFTBDY 

735 


340 

CONTINUE 

AFTBDY 

736 



CALL INFOPLTCOf IMAX.XnUM»l»YOUM,l,XMIN,XMAX,YMINfYMAX» 

AFTBOY 

737 


1 

1.0fOfHEOfOfHEO>0>SIZElfSI7E2*O.Of2.0) 

AFTBOY 

' ‘ 738 


350 

CONTINUE 

AFTBOY 

739 




AFTBDY 

' 740 


C** 

XI - LINES 

AFTBOY 

741 


C** 


AFTBDY 

742 



KK - 0 

AFTBDY 

743 



00 400 I « IflHAX 

AFTBOY 

744 



IF(I.EQ.IHAX) KK * 1 

AFTBOY 

745 



00 38C J ■ IfJHAX ... 

AFTBOY 

746 



XDUH(J) ■ XdfJ.ISTN) 

AFTBDY 

747 



YDUM(J) • YdfJ.ISTN) 

^AFTBDY 

748 


JIQN RUN 


CARDS ENCnUNTEREO IN INPUT UPDATE 1V2-V52. 
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380 CONTINUE 

CALL lNEOPLT(KK#JMAX#XDUMf 1»YDUM#1.XNIN#XNAX#YHIN»YMAX# 
I 1.0i0#HE0»0#HE0#0#SlIEl#SI7E2»0.O#2.O) 

^00 CONTINUE 
C^ + 

A80 CONTINUE 
500 CONTINUE 

C^^ ALL FUNCTIONS COMPLETE 

----- 

900 CONTINUE 

RETURN ' 

1000 F0RMAT<16I5) 

1020 FQRMAT(7HIHAX « #I3>10H JMAX « #T3»10H KHAX - #13) 

1040 FQRHAT(2I5»F10.0) 

1060 F0RHATU7H~ Y-2 PLOT I ■ #I3#7H X - #E11.5) 

1080 FCRHAT(17H~ X-Y PLOT K - #13. 7H Z • #E11.5) 

C + + 

END 


AFTBDY 749 
AFTBOY 750 
AFTBDY 751 
AFTBDY 752 
AFTBDY 753 
AFTBDY 754 
AFTBDY 755 
AFTBDY 756 
AFTBDY 757 
AFTBDY 758 
AFTBDY 759 
AFTBDY 760 
AFTBDY 761 
AFTBDY 762 
AFTBDY 763 
AFTBDY 764 
AFTBDY 765 
AFTBDY 766 
AFTBDY 767 


CORRECTION IDENTS ARE LISTED IN CHRONOLOGICAL ORDER OF INSERTION 

AFTBDY 


S ARE LISTED IN THE ORDER OF THEIR OCCURRENCE ON A NEW PROGRAM LIBRARY IF ONE IS CREATED BY THIS UPDATE 



YANK$$$ AFTBDY 


B. Job Control Cards 


A list of job control cards needed to run this program are 
included in the comment section of the main progrcim AFTBDY. 

C. Input Explanation 

The instructions for input preparation are included in the 
comment section of the main program COORDC. 

D. Sample Input 

Sample input listing follows on the next page. 

E. Sample Output 

Sample output follows sample input. 
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^ < 

O CL 


Q 

QC 

o 


CQ 

o 


UJ 

H 

Z 

UJ 

U 


:S 

o 

Urn 


X 

u 

cei 

< 


z 


o 

o 

u 


to 

UJ 


Q 


Z 

< 


z 

< 

cc 

H 

O 


o 

z 


z 

o 


! — 

u 


> 

CQ 


z 

O 


:S 

^ Q 

UJ co 

O Q > 
ci O ^ 

Q. U Q 



I 
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NOTE; WRITE NUMBERS 10, LETTERS I ^ \J G i C. SYMBOLS 




LANGLEY RESEARCH CENTER 

FORTRAN - DATA CODING FORM 


PROGRAM NO 


CODED BY 

DIVISION 


-SECTION- 


JOB ORDER. 


STATEMCNT 

NUMBCri 


-CONTINUATION 


FORTRAN STATEMENT 


12 3 4 5 


I I.. I I- 


{ J .! ...I - 


l -U .1 .1 


8:2 . 8 


8;6^ 4 

9:0 i O: 


„ 1 -_I 


.1. _J L — 1. . -1- 1 1-. — J— -I — i i. __L — u -i i - — 1 ~~ i — i — I — \ l. I - - i 


J 1—4 


l— J_P_J_.*-1 -Pj I I— l__i I L—X L. .1 — i_ 

8.I.3.I. 


7 4 111, 

L .. L 1 -. 

99135, 

1 I 1. .-J L 

5 1 2 0 0, 

I I 1 l-,_l L L__L_ 

3.5A6.8, 


_1 I L 1 L 1 U. 


-4 X — I J 1 1. -.-4* - 4 1 


-U.4 .i-4- 


1 1 1_—4— I 1 1 1 I L_ 


_l 1 J 4 1 


J, 1 L— li_l I I I J L J I I L 

i._ 1 L I i..a I 1 1 1 1 

1 1. _.L. , .1. -.-L- L- . -_l 1 L L 


18028, I 

23416 , , 

75200, I 

7 5 3 1 4, 

i -.1. -1. . 1--1- L 

2 5 0 4 8, I 

2,5A4,8, ^ I ^ 

JAJA~Aa x,-,!.,.... , I . . . 

0, 0| I J I I I I . 1 ...J 

2, 3, 4, 1,6! , , ., ...1. , 

1.8.0. 2.8| , , ^ I , , ^ 

3.8. 4,5, 3| , , , _, ! 

3, 5, 4, 6, 8| , .,_L .1 .1. 

5. 1.2. 0. 0, , . 


^ L 1 J 1 i_ 1 — I 1 I I J _ 


J 1 1 -.1 J ..._1 t i_ 




__L.,4__i 




1... 1 .. • . I...J l..„ J X i_._l L I ...J L..._L... 

‘ J .! _ I - I J ! 1 i L_ 


L..4 _ .1 - -1 




--_1- - } _.t J. ! 


,_.i I I i i ! !.. 


L 




-J I L„JI J„... 
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NOTE; WRITE NUMBERS 10, LETTERS 1 p' U G ^ C, SYMBOLS 




PROGRAM NO 

CODED BY 

DIVISION- 


LANGLEY RESEARCH CENTER 

FORTRAN - DATA CODING FORM 


PAGE. 


SECTION- 


JOB ORDER- 


GTATXMCNT 

NUMUILU 


-CONTINUATION 


FORTRAN STATEMENT 


2 3 4 .5 6 7 8 9 10 11 .12 13 14 U }6 17 10 19 20 21 22 23 24 25 2§ 27 28 22 3Q 31.32 33 .34. 35 36 37 .38 .32 40 41, 42 43 44 45 .46. 4Z_^a_49 J0_^i_52_33 J4-3i.^6-53L3fi_52 60 .61-62_63-4 


. L -1- 4- 


x_u-L- 


__4 — I — I — i \ 


.!J .?• .7i.A-3l .2] L_„J L L — I 1 i- — J— -I 1—4 1 — I I 1 — 4--.1 — 1 i 1 — U- -1— l . 

ll. 3l . ti.. S _.8l 1 l,.7J L. I L I 1_._ 1- -J— 1 J J u— 4 -I 1 - -i 1 -I 1 1 \ 4 4 4—... l- 

,..li . Qi Xjj-Qi—QL. Qj 1 1 J I — 1 — 1—1 — I — 1 — L — I — I — t — u — J — I — i — I — i — 1 — I — I — i. — i,. 

1 _..L L 1 J I J I 1 1 1 1 1 1 4 J 1 4 1 1 1- 

I — 1 1 I 1.-... I 1 I L.Qj . IJ. 3 L-jL-la.-A i_ 6 1 _3j 1 1 I I 1 1 J 1 I I J I a I _ 

I 

. 1.. 1 L-i- 1 1 iJL-J4._^L34..34.AuQi-J 31 i 1 J L L 1 1 U-4 1 1 1 1 4- 

I 

.1. .1 - i _.; J I L I 1 1 3[ I I I I I , I 1 1 I J J 

I 

4 I I.. ! l.. _4 i 1 1 ...4 1 I J J J 1 J J J I i I- -I J_ 

... i 1. .. 6. 6. 4. 0. 6. 3 . . 




1 1 i 1 i » 1 1 i- I- --4 1 — ..4 1- -4 - 


J— a— 1— 


1 1 1 


_. 1 L-. 1 1 . 




I J. ...L 1 1 i ?L.5i Ai Pj Oj J . L..J I 1 J.. 

L....L-_.. 2, 5j.2, .4, 4^6^ 

, 3 . 0, 1 4 4 0 0, , 

I 1 1 L. I 1— i-.—l J L—l 1 4 I I i..„l i— i 1 1— 


-.1 


-i. -L- 


.t— i * I 4. 7. 4 . 6. 3 [ . ..J . . .. I — . — . J . I . 

I 4i _?_i 1 l ?j _5i_Qj_0l.0l I 1 L 1.1 — i — 1 I 1 1 I 

- ^ l\ .'^L .1. ...1__I._1 _ 1 __1_ ... 

. 1 .5. , .•.,.4| .1, 4^ 4^._Pi _ ...... 1 

A,^,.h. 2i .6,.. 3i ......i._..i... j._. 1.. 1. .1...... I 

! .^1. .* I. .^.J. .Pi Pi I. - ...J_ ._J J...I. .. 1 ! I...I .1. .1.. 


J .. . ! I - . i 


I , 7, 6| 6, 4, 0, 6, 3| 

I , , 8, ., 4| 8, 6, 4, 0, Pi , , 

! 1 ,9, .,3| 3, 7i 4, 6, 3; , , 

I ,1,0, .,2|1, 1, 4, O^Oi , , 

I. .........i 

I ,1,2, .,0|0,0,0,0,0| , , 


L...J L.1. 

, I . 




_J L 1 1 1... -1 _l _t. .. I. _ i ..i . » 


,;;A-L-in;.':io.y F..nn Oj (M/,R G'j) 


NOTE: WRITE NUMBERS 10, LETTERS I p' U G ^ C, SYMBOLS / . , 
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-f 


PROGI^AM NO 

CODED BY 

DIVISION 


LANGLEY RESEARCH CENTER 

FORTRAN - DATA CODING FORM 


-SECTION. 


D 

P, 

JOB ORDER-. 


STArrMCNT 

NIIMUER 


12 3 4 5 

. 1 


.! .... !. .. L . .1 . 


1 > ' 

I, . ! . I I 

} . .. L I- -I. t 

2 


.! L. - 1. .1 

I I I I 

r 

I I I .1 i 

I .1.1 1.6 

l .. .. l . . ! .1 ^1 

I .. .L . . I . 1 

I I. .1. ». 


I I I I 


W' 


-CONTINUATION 


FORTRAN STATEMCNT 


7 , 8 9 10 11 12 13 14 IS 16 17 18 19 20 21 22 23 24 25. 26 27 28 29 30 31. 32 33 34. 




- ! 


. „ 1 . . j ... 

1 1 1 1 1 J_ J 4-.. 1 1 1- -1 _-4- 1 1 1 111 

• 1 1 1 

1 4 - 1 - 

4 L._ 

4 1— 


J 

1 

1 

1 


1 

1 

1 

J 
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4 
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. 4 
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1 1 i 1 

111 
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_4 


4 

_a 


1 


1 

L_ 

1 1 

1 

1 


.2^.21-,. 

i. _1 —I. J.. ...L. J 1. 4 ._l 1 1 -1 1 1 1 

. f 1 

1 _ 1 __ 1 


._A --J- 

o_ 





1 

_0._ 




._L_ 


1 


1 . 

2,1„ 

i 

1 1 1 j 1 4 4, •! 0| 1 , , , 1 . 1 1 _ J. I I 1 11 

1 1 1 ■ 

-_J i__i_ 

._4 4_ 

1 t 




1 
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1 . 

_L_ 




1 i 

1 

1 

1 . 

.. 1 1... 

1 

1 1 1 J 1 1 • 1 1 1 1 1 1 1 11 1 1 1 L I 

1 1 t f . 


t_ . 4 





1 
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1 

1 1 

1, i ... 

1 

1 1 1 ! 1 1 ^ 1 • 1 ® 1 1 1 4 1 . L 1 1 . 4 , 4 1 1 1.1. 
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_i 1_ 

U- 







1 
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JL 
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l_ 



2,. li .. 
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2 2i 
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I 
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1 L 
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J L- 
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- .10051E + C2 .35475E + 01. 0. .10031E+02 .41250E+0l._, . .0. .10051E + 02 .47475E*0l- 

. .10051E+02 .54144E+.01 d. a0Q51E+02. .61245E+01 0. ,1005XE*02 .68754E*0l 

. .10051E + 02 .76641E + 01 tt. ^10051Et02 .84864E+01 0. .10051E + 02 ..93373E+01 

. ^10051Et02 .lQ211Et02 Q. .10051Ef02 .lllOlE+02 0. .10051 £+02 --.12000E+02 


I • 1 J ■ 26 (X»Y>2» .. 

0. .1500CE+02. _..0. . . 0. .15000E + 02 .30146E+00 _0. .15000E + 02 .61140E+00 

0. .15000E + 02 .93746E + 00 ..0, ..15000E + 02 .12864E + 01 0. ,15000E+02 *16641E+01 

0. .15000E+02 ...20754E+01 0. .15000E+02 .25245E+01 0. .15000E+02 ..30144E+01, 

0,. ...15000E + C2 . .35475E+0L, 0. .1500QE+02 .41250E+01, 0. .15000E + 02 .47475E+01 

0. ..15000E + 02 .54144E + 01 .0. .15000E+02 .61245E+01.... .0. ..15000E + 02 .68754E+01- 

0. ,.,15000E + 02 . .76641E+01 _0. .15000E + 02 .848648 + 01 0. ..15000E + 02 .93375E+01_ 

0._ . .15000E + 02 .10211E + 02 0. .15000E + 02 .lllOlE+02— 0. .15000E+02 .12000E+02 


I ■ 6 J ■. 1 (XiY.2) 

.. 18000 E +02 . 123878+01 

. ieOOOE+02 . 123878+01 

. 180008+02 . 123878+01 

. 180008+02 . 125878+01 

. . 160008+02 . 125878+01 

,. 180008+02 . . 123878+01 

. . 180008+02 . . 125678+01 


0. .180008+02 .125878 + 01 .301468+00 .180008+02 .125878+01 .61140E+00 

.937468+00 .180008+02 .125878+01 .128648+01 .180008+02 .125878+01 — .166418+01, 

.207548 + 01 .180008 + 02 .125878 + 01, . .252458+01 .180008+02 .125878+01, .30144E+01 

.354758+01 _ .180008+02 .125878 + 01 .412508 + 01 .180008 + 02 .125878 + 01 .474758+01 

.541448+01 .180008 + 02 .125878+01, .612458+01 .160008+02 .125878+01 .687548+01- 

.766418+01 .180008 + 02 .125878+01 .848648+01 - .180008+02- .125878+01 .933758+01- 

.102118 + 02 .180008 + 02 .125878 + 01 .111018+02 . 180008+02 ___.125878 + 01 .120008+02 


1 - 6 J • 

. 160008+02 
. 180008+02 
. 180008+02 
. 160008+02 
. 180008+02 
, . 180008+02 
. 180008+02 


,6 <X,Y.Z) 

. 222178+01 
. 222178+01 
. 222178+01 
. 222178+01 
. 222178+01 
. 222178+01 
. 222178+01 


0 . 

.937468+00 

.207548+01 

.354758+01 

.541448+01 

.766418+01 

.102118+02 


.180008 + 02 .222178+01 _ .301468+00 .180008+02 .222178+01 .611408+00. 

.180008+02 . .222178+01 .128648+01 .180008+02 - -.222178+01 .166418+01 

.180008 + 02 - .222178+01 .252458 + 01 .180008+02 .222178 + 01 .30144E+01 

.180008+02 .222178+01 .412508+01 .180008+02 .222178+01 . .474758+01 

.180008 + 02 .222178 + 01 .612458 + 01 .180008+02 .222178 + 01 .687548*01 

.180008 + 02 .222178+01 .848648 + 01 .180008+02— .222178+01 ..933758*01 

.180008+02 . .222178+01 .111018+02,, .180008+02 .222178+01 .120008+02 


1- 6 J, -.11 (X>Y#Z) . .. . - - 

_. 16 0008 + 02 .389268 + 01 .... 0. .180008 + 02 .389268 + 01 .301468*00 .180008 + 02.^ .389268 + 01 .611-408+00 

-.180008 + 02 .389268+01 .937468 + 00 .160008*02 .389268 + 01— .128648 + 01 .180008 + 02 .389268+01, -—,.166418*01 

.180008+02 . .389268 + 01— .207548+01 .160008+02 .389268+01 .252458*01 .180008*02 .389268+01 .301448*01 

... 180008*02. _ .389268+01 . 354758+01— _.l 80008+02 .389268+01^ .412508 + 01 .180006 + 02 .389268*01 .474758+01 

.180008 + 02 .389266 + 01 ... .541448 + 01 .180008 + 02— .389268 + 01. .612458+01 .180008+02 -. .389268 + 01— .687548+01-. 

.180006+02 .389266+01 .766418+01 .160008+02 .389268+01 .848648+01 .180008+02 .389266+01 .933758+01 


" WK5T ■: ;* ' iwi r ^ * srr/T^ •^r*.’? 


18000E + 0Z •38926E + 01 .10211E*02 -.-UB000E>02 -♦38926E*0X •11101E*02 .IBOCOE^OZ •3e926E*0l •12000E^02 


I 6_ J _-.^6 iXsiMll 

.leDOOE^OZ . . .66^09E+01 . ^0, . .... .ie000EtO2, «664Q9EfO] «30146Ei00 -..•IBOOOEf 02 466409E401-. «611A0E«Q0 

.lBOOOE+02 .66409Lf01 ^937A6EJ^Qi) ♦1P000E+02-. •B64Q9Et01 .a2a64E.^01 •1B000E^03 .66409E401 *166^1E>01 

.16000E + 02 .66409E + 01 .2075AE + 01_ .18000Ef02 .66^09E + 01 .252A5E+01 ..lBOOOE+02 •66409E+01 . •30144E+01 - 

.18000E+02 .66^09E^C1. . • 35475E+01 . _ .18000E*02. .66^09E + 01 •41250E+01 •lBOOOE+02 — .66A09E^01 .A7475E>01 • 

.18000E + 02 .66^09E + 01 .5^14^EJ:i)l ^IBOOOE-fOZ. *664Q9E + 01 .61245E-i0i .lflOOOE + 02 .66409E«01 •6B754E>01 

.I6U00E^02 .66^09E401 •76641EfQJ .1B000E>Q2 .66409E+0I .B4864Ef 01__.1BOCOE + 02 •66A09E + 01 •93375E>01 

.18000E+02 _..66^09E+01 •10211Ef02 •16000E+02 .^66409Ef01 .11101E*02 .IflOOOE+Oa .66^09E+01 ♦12000Ei02- 


■ 6 J .-. .21. (XfYiZ) 

.18000E + 02 . .lC^66E + 02 . 0. „ 1 POOQE+02 .10466E+02 .. .30146E+00 .IBOOOEfOZ .10466E+02 .61l^0Ef00 .. 

.1800UE+02 .10^66E+02_ •93746Ef00 . .. .leOOOEtOZ *10466E+02 . .. .1286^Ef01 .lBOOOE+02 .10^66Ef02 16641E-♦01-. 

_ .18000E+02 .10466E + 02 .- .2075^£ + Ql .18000Ef02 .10466E + 02 .25245Ef01 .1BOOOE+Q2 •10466Ef02 .30144E-f01- 

.10OOOE+O2 _..10466E+02 .3b475E*01 _.18000Ef02 .10466Ei0i .41250Ef0l .18000E+02 .10466Ef02 . ...-.•47475Ef 01... 

.ieOOOE + 02 a0^66E + 02. .5414 4E+01 .16000E + 02 .10466E + 02 . . .61245E+01 .18000£f02 .10466E>02. .. .68754E+0L.- 

.18000E + 02 .10466E+02 _ .. .76641E + 01 .lP000Ef02 .10466Ef02 .a4864EfOl .IBOOOE^OZ . - .10466E+02 -.93375E401 

.. . ..18000E + 02 .10466E + 02 .10211E+02 . .lP000Ef02 ... .10466E + 02 .lllOlEfOZ .IBOOOE^OZ .10466E+02 .12000E402- 


.. I • 6 J - 26 (XfY#Z) , 

.laOOOE^-OZ .15000E+02 . 0. .ieOOOE-^02 .15000E^02 . ...>.30146E+00. . .lBC00Ef02 .l5000Ef02 - .61140E+00.. 

.18000E+02 .1500CE+02 .93746E+00 .lPOOOE+02 .15000E + 02 ... .12664EfOl .IBOOOE^OZ -. .15000E+02 .16641E401- 

.ieOOOE+02 .15000E + 02 .20754E + 01 .. .leOOOEtOZ .15000E^02 .. .25245E+01 .IBOOOEiOZ .15000E+02 «30144E401 

.18000E+02 .15000E+02 .3b475E+01 .lPOOOE+02 .15000E^02 .41250E+01 .1B000E>02 .15000E + 02 - - - .47475E*0l - 

.18000E+02 .15000E+02 .54144E+01 .18000E+02 .15000E + 02 .61245Ei01 .16000E+02 .15000E+02 .. .6B754E+01-. 

.16000E+02 .. ..15000E+02 . .7664lEf01 .lPOOOE + 02 .15000E + 02 ..a4B64EfOI .lB000EfO2 .15000E + 02 *93375Et01 

.18000E + 02 _ ..15000E+02 .10211E + 02 .18000E+02 .....15000Ef02 . .. .lllOlEiOZ .IBOOOE^-OZ .15000E+02 - - .12000E402 > 


I • 11 J ■ 

.36000E + 02 

.36000E+02 

»3tOutE'*^o2 

.3t00uE+02 

,36000E+02. 

»3600uE4t2 

.36000E+02 


1 (X>Y>Z) 

•25174E+01 
.25174E+01 
•25174E+01 
.2!?i74E + Ul 
• 25174E + 01 
.25174E+01 
•25174E+01 


0. .36000E+02 

.93746E + 00 .3^00^F♦n2.. 

.2075^E^0l .36000E*02. 

• 3547*^F + 01 .36OO0E + O2 

. .54144E+01 .3EOOOE+02. 

• 76ft41lP+01 .3ftOO0E + O2_ 

.102llE+n2 __ .3f000E+02. 


•25174E+01 

.25174E+01 

•25174E+01 

.25174E+0L 

^25174Ei01_ 

.25I74E+01 

• 2517.4Ei01. 


.30146E + 00 

_.1?864E+01 

.25245Ei01 . ... 
.j41250E + 01_. 

..61245E101 

..B4864E+01 

.lllOlE+02 


•36000E^02 
.36C00Ei02 
•36000E402 
..J6000Ef02. 

.•36000E>.02. 

„.36Q00Et02. 

•36000Ef02 .25174E401 


• 25174E + 01-.- - .61140E+00 

• 25174E401 ♦16641E+01— 1 

.25174E^01 . .. .30144E+01. ^ 

.25174E + Q1 .4.7475£ii)] 

.25174EtOJ .6a754Ef01 


„.25174E + 01_ 


.^93375Efat.. 

•12000E>02 


11 J ■ . 6 (X,Y>Z) 

.36000E + 02 ...33921E + 01 .0... . .36Q00E:f02 .33521E10] .*30146Ei0a .36000E+02 .33921Ef01., — .£1140E^0Q 

.36000E^02 .33921E + 01 , . .93746E^00 .36000Ef02 .33921E+01 .12Q64E + 01 .36000E+02 -.-.33921E + 01 .16641Ei01 

.36000E^02 ... .33921E401 . .2075AEf01 .3600nE + 02 .33921E+01 .25245E*01 .36000Ef02 .33921E+01 .30144E401 

.36000E+02 .3 3921E401 .35^75E+ai .36000E*Q2 .33921Ef01 .41250Et01 .36000Ef02 .33921E + 01 .47475Ei01 

.360008^02 .33921E«^01 .54144Ei01 *36000E^02 .3392164^01 .61248E401 .36000E>02 .33921E>01 .68754E401 



•36000E + 0L_ ..33921E + 01 .76641E + 01 ..36000E+02 .33921E+01 .8486AE + 01 .36000E«02 - . .. .33921E^01 .93375E+Q1 

.36000Et02 •33921Et01 .10211EtQ2 ^ED0DEi02 4 33921E10L •.UlAl£,t.02 .36000E*02.^ »33921E«0] «-12000E^02 


J..! ll_.tXiY>2J. 

.36000E + 02. .A910CE + 01 0. *36Q00E*02 .A9100E»0L ,30146E*00 .36000Et02. -.49100E + 01 .61140E+00 

.36000E+02 ,<i9100E+01 .937'i6E*00 .36000E+02 ^(49100E*01 .12864E + 01 .36000E+02 - .A9100E+01 .... .16641E+01 

.36000E + 02 ,^9100E + 01 ,2075%Et0J »3EOOOE*02 .49100Ef01 .25245E+01. .36000Et02 .49100E+01 ^3014^E+01 - 

.36000E+02 ,'i9100E+01 .35475E+Q1 .3EO00E+02 .A9100EtOl .412S0E + 01 .36000E+02 .49100E + 01 .47A75E+01 

.36000E + 02 .49100E+01 .54144E+01 .3t000Et02 »49100Et01 .61245E+01 .36000E+02 .49100E + 01 .68754E+01 

.36000Et02 .49100E+Q1 .76641E + 01 .3t0O0Et02 ^49100E*01 «84a64E + 01 .36000E + 02 .49100E+01 .93375E+01 

.36000E+02 .49100Et01 .10211E + 02 .360Q0E±02 ..49100E+01 .11101Et02 *36000E*02 .49100E*01 .12000E+02 


.36000E+02 __ .74066E+01 ... 0. .36000E+02 .74066E*01 .30146E + 00 .36000E+02 .74066E+01 -..61140E400 

.36000L+02 .74066E + 01 .93746E+00 .36000E+02 .740666*01 .12864E+01 .36000E+02 — .-.74066E + 01 .166416*01 

_.3600C£*02 .740666*01 . . .20754E+01 .360006*02 .740666*01 _ . .252456+01 __ .360006*02 .740666*01 .301446*01 

.360006*02 .. . .74066E+01 .35475E+01 .360006*02 .740666*01 .412506*01 .360006*02 .740666*01-. _.47475E*01 

.360006+02 .. .740666*01 ...541446*01 .360006*02 .740666*01 .612456*01 .360006*02 .740666*01 .687546*01 

.3600UE+U2 .. .740666*01 .766416*01 .360006*02 .740666*01 .848646*01 .360006*02 .74066E+01 .933756*01 

.360006*02 . -.740666*01 .102116 + 02. .360006*02 .740666+01 .111016*02 .360006*02- -.740666*01 1-.12000E+02 


.360006*02 
- .360006*02 
. .36000E+02 
.36000E+02 
.36000E+02 
. .36000E+02 
__ .36000E+02 


-21 (X,Y>2) 

.10B82E+02 
.10882E+02 
.108B2E+02 
.108826*02 
.10882E+02 
.108826*02 
.108826*02 


0. .360006*02 .108826*02 .. .301466*00 -.360006*02 .108826*02 -.611406*00 

.937466*00 . .360006 + 02 .108826+02 . ...128646*01 .360006*02 .108826*02 .166416*01 

.207546*01 .360006*02. .108826*02 .252456*01--.. .360006*02 .108826*02 -.301446*01 

.354756*01 .. .360006*02 .. .108826*02 . .412506*01 .360006*02— .108826*02 .474756*01 

.541446 + 01 . . .360006*02 . ... .108826+02, .612456*01 .360006*02 .108826 + 02 .687546*01 

.766416 + 01. ..360006*02- .108826+02 ..848646 + 01 - .. ..36000E+02 .108826+02 -.933756*01 ; 

. .1021 16*02 .360006*02 .108826*02 .111016*02 .360006+02 .108826 + 02 .120006*02 


I » 11 J « 26 (X>Y,2) . 

.36000£*02 .15000E + C2 0. .360006*02 .150006*02 .301466*00..- .36CC0E+02 .150006*02 .611406*00 

.3tojUE+u2 .150016*02 .937466*01 .360006*02 . .150006*02 .120646*01 .360006*02- .150006*02 .166416*01 

,3t0o0E*O2 .150006*02 .207546*01 .360006*02 _ ,.150006*02 -..252456*01 .360006*02 .150006*02 .301446*01 

. .360O0L+02 .1500O6 + 02 .354756*01. , .360006*02 .150006*02 .412506 + 01— .360006+02 .150006 + 02 ^.74756+01 

.3600oE+02 .15000E + 02 .541446 + 01 .360006+02 - .150006 + 02 . .612456+01 . .360006+02 .150006+02 _ .687546+01 

, .360006*02 .150006*02 .766416*01 . . .360006*02. .150006*02 .848646*01 .360006*02 .150006*02 .933756*01 

.360006+02 .150006*02 .. 102116*02 .360006 + 02 .150006*02 .111016+02- .360006*02 — .150006 + 02 .120006 + 02 


I ■ 16 J 

.54000E+02 
.540006*02 
. ,5400uE*02 
.540006*02 


.377606*01 0. .540006*02- .377606+01 ,,.301466*00 .540006+02 .377606+01- _.61140E+00_ 

.377606 + 01 ..937466 + 00-. .540006+02 .377636+01 .128646+01 .540006+02 ,377606 + 01 .16641E+Q1-. 

, .377606 + 01 .207546*01 .540006*02 .37.7606 + 01 252456+01 .540006+02 .377606*01 ,30144E-*01__ 

,377606+01 .354756+01 .540006+02 .377606+01 .412506+01 .540006+02 .377606+01 .474756+01 


.t»^000E+02, - .37760E + 01 .S^HAE+Ol •S^OOOEtOa .37760Ef01,_ „_.61245Ef01 •5^000Et02 •37760E+01 -•6875AE+01, 

• 5^OOGEi02 •37760E+01 ..76641E+QJ .5f0Q0E*02. •37760Ef0L ^ _^84B64Ef0l •54000Ef02 .37760E>01^- <-.93375E+01 

_«5^000£t02._ .37760E+01 _.10211E + CI2 •5^000E±D2 •3I760Et01 ♦1110lEifl2 •5^0C0Ef02 •37760E*01 UaOOOE^O^ 


1..- . 16 J 6^ (X#Yi2) 

• 5^000E + 02 .^5626E + 01 0. ... . 54000E^02 . ..45626Ef01 .30146E^00 •54000E+02 .. . •^5626E+01 . •611A0E+00 

_.54000E+02 .45626E+01 .93746Ef00 .54000Et02 .45626Ef01 .12a6^Et01 .54000E+02 .A5626E + 01— .16641E401 

._.«5<»000E+02 .<t5626E401 .20754Ef01 .5^000Ef02 .^5626Ef01 _*25245E + 01 •5^000Ei02 ...... .445626E+0] .301A4E+01. 

.5^000E + 02... .^5626E*01 .35475Ef01 *5 4000Ef02 •45626E401 •41250E + 01 •54000E+02 .4 5626Ei01 •47475Et01- 

.54000E+02 ..45626E + 01 .64144E + 01 .54000E + 02 •45626E*01 •61245E+01 .54000E+02 ,45626E+01 .68754E+01 

_.54000E + 02 .. . .45626E+01-.._ .76641E + 01 .5 4000E + 02 •45626E^01 •84064EfOl ^^_.54000E+02 •45626E+01 .93375E+01 

..64000E+02 •45626E+01 .. ..10211E + 02 .54000E + 02 a45626Ef01 ..lllOlE+02 .... .54000Ef02 .45626E + 01 .12000E+02 


I • 16 J ■ 

.54000E + 02 .. 

• 54000Et02 ... 

. .54000Et02 

. •54000E+02 

•54000E+02 
. «54000Ef02 
.. •54000E + 02 ; 


11 (X#Y#2) 

.59276E401... 

. •59275E401... 

•5927&E401 

•59275E401 

.59275E+01.. 

•59275E401 

.59275E+01 


0. .54000E402 .59275E401 .30146E400-- •540C0E402 

..93746E400 . . ^ 54000E402 .59275E401 .12864E + 01 •54000E402 .. . 

.20754E + 01 ... .54000E+02-... .59275E401 . « .25245E401 , -•54000E402 . 

• 35475E + 01 - •54000E402 . . .59275E401 .41250E401 •54000E402 ... 

.54144E401 .54000E402 .59275E401 .61245E401 . •54000E402 .. 

• 76641E401 .54000E402-. .59275E401 .84864E401 .54000E402 

• 10211E402 ......54000E402., . •59275E401 •11101E402 •54000E402-- 


• 59275E + 01 ♦61140E400 

..59275E401 .16641E401 

•59275E+01 •30144E401 

.•59275E401 - . ♦47475E401 

• 59275E401 -.68754E401 

.59275E+01 •93375E401 

• 59275E+01 a2000E402- 


1-16 J • 16 (X#Y#2) 

•34000E402 .ei722E401 0« 

•54000E+02 .81722E401 .93746E400 

•54000E402 .81722E+01 .20754E401 

.54000E+02 ..01722E4O1 •35475E + 01 

. .54000E402 ... .81722E401 .54144E401.... 

t:<«^UUUt4u2 . •&1722c4Ci. •766^1f4^l 

.!/40vCb + t*^ .ftl72^L + Cl . + 


. «54000E402 •81722E401 «30146E400 

• 54000E402-. •81722E401 •12864E401 

• 54000E402 .81722E401 •25245E401 

•S4000E402 •81722E401 «41250E401 

354C0CE402 ..^ei722Ef01 •61245E401 

.5^0CCE402 .81722E401 .‘J4a64E + Gl 

.f4COCF+02 .ai722H+01 •11101E+C2 


• 54000E402 . .«81722E401 *611408400 

•64000E402 .81722E401 - •16641E40t 

• 54000E402 — .B1722E401 .*30144E401- 

•54000E402 .81722E401 *47475E401 

•54000E402 .81722E401 --.68754E40t 

• 5 4CtCL402 .81722E401 ... . . •93375E401- 

.^4CCoE+02-. -.ei722E + Cl ..-.-*12000E402 


I » lt> J * 21 <X> Y>Z) . . . . 

• P4000L + C2 ..._,11297E + 02 . 0.._ ..*^40008402 ..11297Et02 .3Q146E + 00 .54GCLE4Q2 •11297E4C2 .611408400. 

• 54000E402 .11297E4U2 .93746E400 .540008402.. .112O7E402 .12864E401. .54u00E402 .11297E402 ..... . 16641E401 

.54000E402 .11297E402 .20754E401 .54000E402 .11297E402 .252458401 .54000E402 .11297E402 .301448401 

.54000E + 02 .11297E402 .35475E4Q1 .54000E402 .11297E402 .41250E401 .54000E402 .112978402 --^ .47475E40U 

.54000E402 .11297E402 .54144E401 .54000E402 .11297E402 .61245E401 . .54000E402 11297E402 .68754E401 

...54000E402 . .11297E402 .76641E401 .54000E4.02 .11297E402. ..84864E401 .540008402 -^.112978402 .933758401 

...54000E402 . .. .11297E402 . .....10211E402 .54000E402 .11297E402 - ......11101E402 .54000E402. -.— .112978402 .12000E402 


1 « 16 J ■ . 

. .54000E402 
..54000E+02 _ 
.54000E+02 
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•28752E+01 

.766A1E^01 

0. 

.28752E+01 

•BA86AE401 

0. 

•28752E+01 

. 93375 E +01 

0* 



.28762E + 0i . 

.10211Ef02. 

„-.Q, 

^2B752E+01 ... 

- aiiioiE+oa 

__.o. 

.28752E+01 

.1 7nnnc«n» 

Z - 1 • 

1 

J ■ 

16 (XfYiZ) 








_ 0 . 



.58752E+01 

0. . Z 

0. 

- — ..58752EtoiZ- 

.301A6E400. 

_ 0* 

t5P75?F+01 

f A1 1 4np4nn 

0. 



.58752E+01 

.937A6E+00 

. 0. 

. - . .5B752E+01 

.-.128646+01 

0. „ 

-.58752E + 01 

.16641E+01 

o,... 



.58752E*01,_ 

. .2075AE+01 

...0. 

.5B752Ef01_ 

— •25245E+01 

0. 

-.58752E+01- 


0 • 



.58752Ef01 

•35A75E+01 

0. 

.58752E^01 

•41250E+01 

0. 

.58752E+01 

•47475E+01 


OOOOCOO^; O O O O p O O 1-4 


0. . ,.5B75aE*01 .. .54144E^01 0. .58752E401 .61245Et01 . 0. .5e792Ei01 .68754E*0l_ 

0. ,58752E*01 .76641E*0l 0. .58752E.4Q1 .84864E401 0. .Se752EtO] *9337SE>01 

Ot ,5e752Et0l .10211Et02 0. .58752E*01 .11101Et02 0. *5B752E+01 .12000E^02_ 


• 1 0„5_21_(X*YjZ) 

.10051E+02 0. 0. 

.10051E+02 ,93746Et00 0. 

.10051E+02 ,2075«t0l 0. 

.10051E+02 . .35^75Et01 Q. 

.10051Et02 .5-tl4^E+0l ... 0. 

.10051E+02 .76641E+01 8. 

- -_.10051E + 02 ..10211E*02--.- 0. 


.10051E402 . .30146E+00 0. .10091E«02 . .611A0E400 

.10051Et02 .12B64E*0l. 0* .10091E402 .16641E>01 

.1005lEf02 ._..29245E»0l._. O.— .10091E40i. ._ .30l44E*0l 

.10091E402 .41250E + 01 0. .10091Ef02 .47479E+0L 

.10091E+02 *61245E^01 0. .10091E402 ^68794E+01 - 

a0051Et02 .B4B64E+01 0. .10051E+02 ,93379E*01 

a0051E*02 ail01E402 — 0. a0051E«02 a2000E>02 . 


« 


.1 ' i 


26 (X»Y,2). - 

. .a5000E+02.- .0. 0. a9000E»02 .30146Ei00 Q. 

,._a5000E*02 .93746E + 00 0. ^15000E*0i a2B64E+01 0. 

a5000E+02 . ,20754E*0l ,._0. a9000E*02- .. ..29245E*0L 0. 

. .15000E+02 _ .35475E+01 0. a5000Et02 .41250E+01-- 0 . 

. a5000E + 02 .54144E+0L 0 . a9000E+02 - .61245E + 01 — . 0 . 

.15000E+02 .76641E+01— .0. .15000Et02 .84B64E+01 ._.0. 

a5000E+02 ao21lE+02. _ 0. _ a5000E*02 ^mOlE+02 _ 0. 


a9000E*02 .61140E+00 

.19000E+02 «16641E.t01_ 

.19000E+02 .30144E401 - 

.19000E+02 .47475E401- 

.15000Ef 02 .6879AE401 . 

- .19000E + 02 .93375E+01-- 

- a9000E402 — .J2000E+02- 


I ■ 6 J • 1 (X>Y>Z) 

aeOOOE+02 a2587E*01 0. 

a8000E+02 a2587E+01 .93746E+00 

aeOOOEtOZ a2587E + 01 .20754E+01 . 

a6000E+02 .J12587E+01 .35475E*0l 

. a8000E*02 ..12587E+01 .54144E*01 . 

_ .iBOOOE+02 a2587E*01 .76641E*01 

- a8000E+02 a2587E+01 .10211Et02 


ae0O0E*02 a2587Et01 - - .30146E+00 __ 

a6000E«02 .12587Et01 . a2864E+01 

aeoooE+02 azsezE+oi. .2524 se^oi 

aaoOOEtOZ--. - a2587E+01 -- .41250E+01 

aeoooEtoz azsazE+oi .61245 e+oi _ 

aeOOOE+02 .12587E + 01 .84864E + 01 - 

a8000E + 02 a2587Et01 ailOlEEOZ 


.18000E*02 .12587E401 .61140E+00-- 

aa000E*02 - aZBaZEtOl .16641Et01 

ae000E«O2 .12S87E*01 — .301A4E*01 

-aeOOOEtOZ .129B7E«01 - .47475E+01- 

-aaOOOE*02 .12S87E«01 - ..68754E+01- 

-.16000E+02 a2S87E^01 .93375E*01 

_.18000E*02 .129B7E«01 .12000E«02- 


I • 6 J - 6 (X»Y»Z» 

aeo00£+02 _ .22217E+01 

aaOOOE+02 .226A9E>01 

aBOOOEtOZ . .22849E«01 

aBoOOE+02 .22945E+Oi 

aeOOOE+02 .22982E + 01 

. aeO 0 oE*O 2 .22967E+01 

..a8000E+02 .22809E+01 


0. aaOOOEtOZ .ZZAOaEtOl - .301A6E«00 . .ISDOOE^-OZ .2294'2E«01 . .61140E400 

.93746E*00 .ie000E*02 . .22732E+01 ... ,12864E + 01 .18000E+02 _ -.22798E+01- .166410401- 

.20754E+01 .ie000E402 .22890E401 .252450*01 .18000E+02 .22922E40I .301440*01- 

.35475E*01 .ie000E*02 .229630*01 _ .412500*01 . .180000*02 .22975E*01 .474750*01- 

.54144E*01 .1P000E*02 .22983E*01 . _ .61245E*01 . .180000*02 .22979E*01 .68754E*0l- 

.76641E*01 .180000*02 .229420*01 .848640*01 .180000*02 .228960*01 .933750*01 

.1021lE*02—._. 180000*02 .226290*01 .111010*02 .180000*02 .222170*01 .120000*02.. 


I - 6 J - 11 <X»YiZI - 

.180000*02 .389260*01 . 0. __ .180000*02 .391810*01 .301460*00 .180000*03 .394080*01 .611400*00 

_-.ieOOoE*02 .396120*01 .93746E*Q0 .160000*02 .397.96E*01___.12a64Et01 .180000*02 .399610*01 .166410*01 

.180000*02 .401070*01 .207540*01 .180000*02 .402350*01 .252450*01 .180000*02 .403420*01 .301440*01 



, .leOOOE+02^ •40^31Et01_ .35475EtOI .IBOOQE.^02 •4Q495EfOX •41250£*01.^: ^lBOOOE+02. •40549E401 •47475E401 

. -leOOOEfO^ .40577Et0J .54144E + Q1 *ie000E+O2 r405.a4E+Ql ♦61245E>i)l aaDOOEm •40566E*01 •60754E*O1_ 

...ieOOCE+02 .40515E+01_ .76641Ef01 .1PQOOE±02^ .^4Q421E+Q1 ^fl4a64E*01 aa000Ei02 •4026lEf01 •93375E401.-.. 

.18Q00Et02__ •40001Et01 .10211 Ej^.Q 2 a£QQ0E±Q2 a95 84E*QJ ail01E4^02 ^I8000E*02 .aE926E*0L a2000E-*02 


I.- 6 J - .16 (X#Y/ZI . .. .. . 

.18000E + 02 ___ .66409E + 0I *16000Et02 ^66632Et01 430146E+0Q .18000E+02 .aEa^SE + Ol .611^0Ei00_> 

_4ie000E + 02 .67060E + 01 •93746Et00 .leOCOE+OE .67249Ef 01 . .IZaa^iE + Ol .18000E+02 .67443E + 01 - .16641E+or^ 

.ieOOOE + 02 .67620E + O1,.__^2O754E+Q1 -.leOOOE.tOZ 467B03E*0J 425245E+01 41B000E+02 .67962E+01 aOl^AEfOl- 

.. .10OOOE+O2 .68103E + 01 • 3ii475E +01_ .. . .. . 1 eO0OEtO2 468219E+01 441250E + 01 ^18000E*02 .6B307E + 01 .47475E+01- 

_,.18000E+02 4 68361E+01 •5414AEt01 .ISOOOEfOZ 468375E+01 .61245E^01 4l8000Ef02 .6e341E*01 •68754E*01 

.16000E+02 .68251E+01 .76641E+01 . .. .. 1 6000E + 02 .68091E>01 .84864E + 01 . .18000E + 02 •67847E + 01 .93375Ei01 — 

.•18000E + 02 .67499E + 01 *1021 IE *02 .ie000E.f02 467025Ef01 .lll01Ef02 .lBOOOE+02 .66409Ef0l ^♦12000E+02_ 


1 • 6 J -....Zl (XjYiZl. 

..ieOOOE + 02 .10466E + 02 ... 0 

. .18000E + 02 . .10504E + 02 . . 

. .1800OE + 02 . .10540E + 02 .. 

.18000E+02 .10573E+02 

.lBOOOE+02 .10592E+02 . 

- •leOOOE+02 .10584E+02 . 

.ieoOOE+02 .10532E+02 


.93746E + 00..,- 
.20754E+01 

•35475E+01 

.54144E + 01 ... 
•76641E+01 
.lOZIlE+02 . 


.•18000E+02 

.ie 000 E 402 

..ieoOOE + 02 

• 1P000E402. 
.18000E + 02... 

• 18000E+02 - 

• 18000E+02 . .. 


..10479E+02 , 
.10516Ef02 
.10552E+02 
-.1O501E+O2 
-•10593E^02 
4l0572Ei02 
•10503E+02 


.30146E+00 .. 
• 12864Ei01_. 
.25245E+01 - . 
•41250Ef01 _ 
.61245E401__. 
.84864Ei01 — 
.tlllOlE^OZ...- 


..18000Ef02.-. 
4l8000Ef02 .. 
..18000E+02 — 
.•18000Ei02-... 
•18000E+02 — 
..18000Et02 ... 
» 180 a 0 E 402 


.•10492£i02. . 
...10526E<*02 
..10563E+02 . 
.•105B8Ef02 
-•10591E+02 - 
•10555E+02 
•10466Ei02 . 


.61140E400. 

-^-.. 16641 E 401 

430144E401 

.. 447475E401- 
-...68754E+01. 

493379E401 

....•12000E402- 


I - 6 J • 

•18000E402 
. .IbODoE^uZ 
.lbOOOE+02 
. •16000Et02 
«lb000E^02 
.ieOOUE+02 
.lbOOOE+02 


26 (X,Y#Z) 

•16000E^02 0. .IPOOOE^OE . .19000E402 

•15000E+02 •93746E+00 •lPOOOE+02. . .15000E402 

.l!>000E+02 •2075^E401 .. .18000E+02 •15000E402 

,.15000E402 . .35475E401_ . . •lPOOOE+02 .15000E402 

•15000E+02 .54144E+01 .1P0O0E402 „ ._.15000E+02 

•15000E+02 .76641E401 .1POOOE402 .15000E+02 

• 15000E + 02 .10211E402 .1P000E402 .ISODOE+OE 


• 30146E+00 . .. .18000E402 415000E402 - •61140E400- 

.12064E4O1 .18000E402 .•15000E+02 > •16641E401-- 

• 25245E401 .. .18000E402 ..15000E402 •30144E401.. 

• 41250E401 . - •18000E402 .15000E402 - . 47475E401. .. 

..61245E401 .18000E402 4 15000E402 •687 54E401^_ 

.84864E401 .18000E402 .15000E402 •93375E401- - 

• 11101E402 •18000E+02 •15000E402 .12000E402 


1 - 11 J • . .1 (X#Y#Z) . ' 

♦ 36000E + 02 •25174E401 0» aeOOOEtOZ *25174E401 . 30146E400 36000E402 425174E4C1 611 40 E 40 

436000E402. ^ 425174E401 . •93746R400 •36000E402 i25174E + 01. ... . 12864E401 _ .. .36000E402 .25174E401. ...16641£4QI_ 

♦ 36000E + 02 •25174E + 01 •20754E401 .36000E402 •25174E+01 •25245E401 _ .36000E402 425174E401 . _ ,..30144E401 

436000E+02 .25174E401 .35475E401 ... •36000E402 .25174E401. . •41250E401 _ .36000E402 .25174E401 447475E401 

♦ 36000E402 •25174E401 •64144R401 .34000E402 •25174E401 .61245E401 . .36000E402 •25174E401 .687546401..,. 

.36000E402 .. .25174E401 .. .76641F401 .36000E402 .25174E401 . _,.84864E401 .36000E402 .25174E401. .93375E401i-. 

.36000E402 . .25174E401 * .10211E402 _ .36Q00E402 .25174E401 .11101E402 .36000E402 .25174E401 .12000E402 


I • 11 J - 

..36000E+02 , 
.36000E+02 


6 (X>YiZ) 

. .33921E401 0. .. . 

•34749E401 .93746E400 


..36000E4Q2- 

•36000E402 


_.34278E401.. 

.34910E401 


..30146E400. 

.12864E401 


.•36000E402.^ 

.36000E402 


__.34544E401_ 

•35038E401 


..61140E400 

•16641E401 


.36000E + 02 ,35l39Et01 .20759E + QI .,36000E*02. .35219E+01 .25245Et01 .36Q00EfO2 .35282E+01 .30144Ef01 

._._.3t.C00£+02 .35329E+01 .35^75E±01 .36000Et02 ^3536'%E.t01 .41250E+01 .36000E+02. .35380E+.O1 .47475E+01 

,36000E*02 _.35^02E + 01 .54144E + 01 ,3fc.00QEt.Q2 i.a3405El;Ql .612V5E.+01 .36000Et02 ^35396E*01 ,6875«E-+01 

.3t0OOE*O2.__...35371EtOl .76641E+01 .36000E+02 .35322E+0J .8^8ME+01 .36000E+02 .35232E+01 .93375E+01 

.36000E + 02 .35060E+01 .10211E.t02 .3E000E±02 .34711E±01 .lllQlE.t02 .36000Ef02 .33921E+01 «12000E+02 


.I...? U J.5-11 (X,Y,Z) 

.36000E + 02 .^910CE.t01 i). .360QOEtO2 .49597E.+01 .30146E+00 .360Q0E+02 .5003'4E>0] .611A0E*00— 

_....3t>00U£+02 .5C<t2^E+01 .937A6EtOQ ,3EO0OEtO2 ^5077«E+D1 .1286^E+0l .36000E+02 .51087E + 01 .16641E*01_ 

.36000E+02 .&136<iEt01 .2075 4Et01 *3E000Et02 .5160^Et01 .25245E+01 .36000E+02 .S1808E+0L .30144E-+01... 

.. . .36000E + 02 .5197gE+0I .35475E+01 .3E000E+02 .52106Et01 .41250E+01 .36000E+02 .52199E + 01 .A7475E+0l._ 

_ . .36000E + 02 ,52253E+01 . . .54144E+01 .3E000E + 02 .52266E + 01 .61245E + 01 .36000E+02 .52231E + 01 .68754E+01 

__,36000E+02 .52136E+01 .?6641E.tO,l ,3E000EtD2 .51956E+0L .84064E + O1 ,36000Et02 .51654E+01 .93375E*01_ 

.36000E+02 .5U63E + 01 .10211EtQ2 .36000E+Q2 .50370E + 01 •lllOlE + 02 .36000E+.02 .49100E+01_1. ...12QQ0E*02,_ 


..A ? 11 J - 16- (XiY#2) 

.36000E+02 . .74066E+01 . 

. .36000E+02 .. . .75322E+01 

. .36000E+02 . .76424E+01 

•36000E+02 - ..77306E+01 . 

• 36000E + 02. . -..77701E + O1 

. .36000E+02 - .77579E+01 

. ...36000E+02 .76176E+01 _ 


0. .36000E + 02 .74508E + 01 .30146E+00 .360C0E+02 .74924E+01 .61140E+0Q. 

.93746E + 00 .36000E + 02 .75705E+01 12864E + 01 .36000E+02 ^76073E + 01 .16641E+01 

• 20754E+01 .36000E+02 .76751E*01 .25245Et01 .36000E+02 .77048E+01 *301A4E*01 

.35475E+01 .36000E+02 .77523E+01 .A1250E + 01 .36000E+02 .77683E + 01 ♦47475E+01 

.5A144E + 01 .3600CE + 02 .77806E+01 .61245E+01 - -.36000E+02 »-77745E+0l ,68754E*01 

• 76641E + 01 ,3fc000E + 02 .77286E+01 .84864E + 01 .36000E+02 .76833E + 01 »93375Et01 

..10211E + 02 .3E000E+02 .75273E+0L .lllOlE+02 .36000E+02 .74066E + 01 .12000E+02 


I - 11 J - 21 <X>Y/Z) - 

.36000E+02 .1O802E + O2 0. .36000E + 02,_ .10907E+02 .30146E + 00 .36000E*02 —.10932E+02 .61140E+00 

.36000E+02 .10956E+02 .93746E + 00 .36000E+02 ..10980E + 02 .12864E+01 . . .36000E+02 , 11004E + 02 .16641E+01 

. .36000E+02 . ..11027E+02 .20754E + 01 .36000E + 02 ,11049E*02 . _.25245E + 01 . _..36000E*02 .11070E+02 .30144E+01 

. .36000E+02 .11089E+02 .35475E + 01 . .36000E+02. .11105E + 02 .41250E+01 .36000E + 02 .11117E + 02 ,.47475E*01 

.36000E + 02 .11124E + 02 .54144E+01 .36000E+02 .. .11126E+02 . .61245E+01 _ .36000E+02 .11122E + 02 ..68754E+01 

.36000E+02 .11109E + 02 .76641E+01 . .36000E+02 .11087E+02 I.84864E + 01 .36000E+02 .11055E+02 .93375E*01 

.36O00E+02 .llOllE + 02 .1021lEt02 , . . , 36000E + 02 .10954E+02 _.11101E + 02 .36000E+02 .10882E+02 ,12000E*02 


I - 11 J - 26 (X,Y,Z) ... . - . 

.36000E + 02 .15000E + 02 0. , 3f000E+02. . ,15000Et02 .30146E+00 „ •36000E*02 •15000Et02__..61140Ef00 

.36000E+02 .15000E + 02 .. .93746E*00 .36000E+02 #15000E+02 . .12064E+O1 .36000E+02 .15000E+02. .16641E+01 

.36000E+02 .1500CE + 02 ,207S^F + 01 ,3f000E+02 . ,15000E*02 .25245E + 01 .36000E*02 .15000E+02 ..30144E+01 

_ .36000E+02 _ .1500CE + 02 .35475F+01 .36000E+02 .15000E+02 41250E*01 .36000E+02 .15000E + 02. .47475E+01 

.36000E + 02 .15000E + 02 .54144E*01._ .36000E + 02 »15000E*02 .61245E+01 .36000E+02 .15000E + 02 .68754E+01— 

_ .36000E+02 .15000E+C2 .76641E+01. . .36000E + 02 .15000E + 02 ..._.84864E + 01 . - .36000E+02 .15000E + 02 .93375E+01-: 

.36000E + 02 .15000E + 02 . 1021 lE + 02 . _ . 3 tOOOE+02 .15000E+02 .11101Et02 .36000E+02 .15000E+02 .12000E*02 


,1 .• 16 j • 1 (X»Y^Z) . 7 .’ '-'' 1 ’ 

,bA000E*02 .37760E+01 0. .54000E+02 .37760E+01 .30146E+00 .540C0E+02 .37760E+01 .61140F*00 



•54000H + 02 ,„.37760E+01 .93746E + 00, .. .S^OODE + OZ ..37760E+01 .1286^E + 01 — ^ .5^000Et02 •37760Et01 «166^1E«01 - . • 

.5A000E + 02 .37760E^01 _ .20754E*01 __.5^000E*02 .37760Ef01 ^25245E>01 .S^OOOEfOZ .37760Ef01 .3014^E>01_ 

„.5^000E + 02 •37760E + U1 _ ^35475E*01 *5A000E+02 ^37760E + 01 ->^41250Ef01 •54000E+02 .37760E401 •^7475E>01 

•5^000Et02 .37760Et01 .S^l^AEfOl .5^000Et02_ .37760E+01 •61245E + 01 •54000E+02 •37760E + 01 .6875^E*0l 

.54000E+02 •37760E + 01 .766^1E + 01 .54000Et02 .37760E + 01 • a486^Et0l .i>^000E+02 — 437760E + 01 •93375E*01 

.5^000E + 02 .3776CiEt01 .10211Eta2 ..54000E+Q2. _.37760Ef01 .lUOlE + 02 ,^54000E+02 .37760E + OX .12000E-*02 


„.5^Q0OE+O2 .45626E+01-.-.0. 5^000E + 02 .^6143Et01 . •30146E + 00 .5^000E+02 .46528E + 01 .6114OE + 0O 

.54000E + 02 .4682^E + 01„. _^937A6E + 00 .5A000E+02 •47059E*01 .12864Ef01 .54000E + 02 .47246E + 01 •16641E^01 

.5^000E^02 .^7395E + 01 .2O754E + 01 .54000Et02 .V751^E>0I .25245Et01 .S4000E + 02 .^^7607E + 0l 430144E+0X 

.5^000£ + 02 .. .47678E + OX . 35^75E + OX , _ . 5 AOOOE + 02 .... •47730E+01 . .^X250Ef0X ... .54000E+02 *47766EiOX 447^75E+0X 

•5<i000Et02 .47787E + 0I .5^14^E + OX _ ...54000Et02 .4779XE.fOX .6X2^5EfOX •5^000E+02 .A7778E + 0X .6875^E40l 

.5^000E+02 . .•A77^XE*OX .7ii6^XE + OX .5A000E + 02 •^7668E + 01 .8486^EfOX *54000E + 02 .47532EfOX .93375E+0X 

*5^000E + 02 .47278EtOX „ •102lXEf02 . f5 A000Ef02._.„.46767Ef0X .XXX0XEf02 .54000E402 •45626E + 0X- .X2000E*02 


I - X6 J •. XX (X#Y#2) . - . 

. •5<i000E + 02 ..39275E*0X . .0. • SAOOOEf 02 .60006E^01 *30X46E + 00 #5^000E+02 •606A3E + 0X-. #6XX40Ef00 

• 5<iC00E+02 •61206E*0X .93746E^09 . .•54000E^02 - . •6X707Ei01 . *X286^E + 0X - #5^000E402 - -•62X54EiOX -- #X664lEfOX 

.5^0001+02 . •625-^8£*0X .2075AE*0X .5^000E + 02 .62889E^0X •25245E + 0X .5^000Ei02 - •63X78E*0X #30XA4E+0X 

.54000E + 02 .oa^XiE+OX •35475E^0l . . • 5 AOOOE + 02 •63599E+01 •^1250E + 0X . 5AOOOE+02 — .6373XE + 0X — • A7^75E^0X * 

• 5^000E + 02 •63809E + 0X .541A^E + 0X ♦ 5A000E + 02 ... .63827Ef 01 •OXEASE + OX •54000E+02 •63777E + 0X ♦6875AE*0X 

• 54000E+02 - ...63641E+01 . .76641E + 0X. . •5^000Ef02 .63387E^01 ; .8^86^E + 0X .54000E^02 •62958E + 0X- .-.93375E401 

.54000E + 02 .62259E + 0X .1O2X1E + 02 .54000E + 02 . .6XX25E*0X .XXXOXE + Oa .5^000Ef02 - •59275E + 0X •X2000E*02 


I - 16 J • 16 CX#y#Z) - 1 — 

.34000E + 02 .8X722E+0X 0. . .SAODOEfOa •8238^E f OX .30X46E + 00-. .54000Et02 .82999E + 0X •6XX40E*00 

.5A00OE<'02 . .U3&81E + 0X . . . .93746E + 00 . 5 AOOOE + 02 .84X36E + 0X .. .1286^E + 0X - .5A000E+02 -- •8466^E*0X •X66AXE+0X 

.54000E + 02 .85X61E + 0X .207.*>AE+01 .5A000E + 02 •8562XE401 , .252A5E+0X ^.5 AOOOE+02 . •86035E*0X .30XAAE+01 ^ : 

• 5AU00E+02 •8639AE + 0X •35A75E+01 •5A000E + 02 .. .86688E+01 .AX250Ei0X .5A000E + 02 .86906E + 0X ,A7A75E^0X 

.. .5A000E + 02 .87039E + 0X .5AX4AE + 0X .8A000E + 02 . •87072E*01 •612A5E+0X _... .5A000E+02 .86989E + 0X .6875AEi01 

.5AGO0E + O2 •86763E+0X .766A1E + 01 .5A000E*02. ..8636XE + 01 .8A86AE + 0X .. .. .5A000E + 02 •85732E + 0X -•93375E^0X - 

• 5A000E + 02 .8A810E + 0X . 1 U 211 E + 02 .54000E+02 •83508E + 01 . .XXXOXE + 02 •5AOOOE*02.-^ .BX722E + 0X —•X2000E+0a 


1 • X6 J 
•5A00UE^02 
.5A00CE+02 
•5AOO0E^O2 
.5AG00E+02 
•5A000E^02 
. .5A000E+02 
.. .5A000E + 02 


• 2X (XiY#2) 

.X1297E*02 0. 

.1XA09E+02 •937A6E+00 

.X1512E*02 .2075AE+01 

.XX599E + 02 ...35475Ef0X. 

.116ABEi02 .3AX4AE+01 

.11627E+02 .76641E+0l_ 

.11A89E+02 .10211E+02 


.5A000E+02 

.54000Et02. 
..5A000E + 02_._ 

•5A000Ef02 

• 54000Et02.^.. 

•54000E+02 

.*54000Et02 


.XX336Ef02 - 

..XlAA 5 Ei 02 
.XX5AAE+02 . 
-•XX62XE + 02 

•X1650Ef02 

..X1597Ef02 -- 
.1XA05E + 02 ... 


•301A6EfOO 

•1266AEtOX 

•252A5E+0X. 

..AX230£t0X. 

• 6X2A3EfOX.. 

«8A88AEfOX.. 

.X1X0XE402. 


-.5A000E+02. 

•5A000E4O2 .. 

. .3A000E+O2 

-.5A000Ef02 . 

..•5A000E+02 

_.5AOOOEt02 .. 
.5AOOOE402.™ - 


-•XX373Ef02 

.XXA79E+02 

-♦XX573Ei02 

.XX637E+02- 

.•XX6AAE+02- 

.•XX552Ef02. 

•XX297E+02- 


^6XlA0E*0a.- 

- .X66A1E+0X — 

•30XAAEA0L— 

.A 7 A 75 E+ 0X-.- 

— .6875AE*0X 

._^93375E40X.^ 
a2000E*02_ 


1 



.54000E*02 _ .X5000E+02 . . 0. . ,5^000E+02 ,15000E*02 - - .301A6E+00 .54000E+02 „ .lSOOOE+02 .6U40E*00 

..5«O00E+y2 .15000E+02 ..937A6E100 .SSOOOE+02 *15000E-t02 .1286^E*fll .54000Et02 »13000E+02 *16641E4Q1..^ 

.54000£*02 .15000E+02 .20754Et01 .5400DE402 .1500QE.t02 .25243E+01 .54000E.+02 ^15000E*02 .3014»E*01 

. .54000E»02,_ .15000E+02 .35475E+01 .5400QE102 .15000E>02 .41250E+0L .64000E*02 .15000Et02 .4747SEf01 

,54000E*02 ,15000Et02 ,54l44Et0l .540OOEf.02 »15000E*Q2 .61245E+01 .54000E+02 .13000E+02— ..68754E+01 

.,.54000£+02 ,15000Ef02 .76641E + 01 .34000Et02 ^15000Et02 .B4864EtDl .54000E*02_: .15000E+02 *93375E401 

..54000E+02 .X5000E+02 .1021JE+02 ,3A000E*Q2 .15000E*02. _4lll01E*02 .54000E+02 .15000E+02 .12000Ei02 


1 • 21 (X>Yj2J 

.72000E + 02 .49315E+Q1 0. .72000E402 .A9314E*01 .30146E100 .72000Et02 .49310E*01 .61140E+00 

.72000E+02 -_.49304E+01 .93746E+00 .72000Et02 ^49294E+01 .12864Et01 .72000E402 .49280E+01 .16641Et01 

.72000E+02 . _ .49260E + 01 .20754E+01 .72000E.+02 449234E»01 425245E*01 472000Et02 449200E+01 43ai44Ef01 

-472000E+02 4 49156E.+ 01 .35475£*01 ,72000E*02 .49100E+01 4 41250E + 01 472000E*02 .49030E+01 <47473E+01 

472000E+02 .. .48945Et01 454144E>01 .72000E+02 448842E+01 461245E*0l 472000E+02 ^48718E*Ol 468754E+01 

. 472000E + 02 448574E+01 .76641Et0l 472000E+02. 4 48406Et01 4 84864E+01 472000E+02 448215E+01 493375E+01 

. .72000E+02 .47999E+01 ,4l0211Et02 .72000E+02 447760E*0l .HlOlE+02 472000E+02 447498E«0l 412000E+02 


I.« ,21 J • . 6 «X>Y»2I . . , 

472000E*02 .56371E+01 0. .72000E+02. 457010E*01 430146E400 ..72000Ef02 .S7481EtO] 461140E«00-. 

472000E+02 .57840E+01 . .93746£*00 472000E+02 458120E + 01 412864E+01 472000E402 . __458337E*01 «16641Et01 

.72000E+02 .58504E + 01 .20754E*01 472000E+02 458625Et01 425245E*01 472000E+02 458707E+01 430l44Et0] 

.72000E+02 , 458753E+01 .35475E+01 .72000E+02 458764E+01 441250E + 01 472000E+02 45B740E*01 .47475E+01 

.72000E+02 . .58682E+01 .54144E + 01 .72000E+02 458586E+01 461245E*01 472000E*02 43B449E+01 .68754E+01 

472000E+02 .58262E+01 .76641F+01 ,72000E+02 .58007E+01 484864E+01 472000E+02._ 457635E+01 493373E+01 

472000E+02 .57136E+01 .10211E+02 . .72000E+02 456290E+01. . 4lll01Et02 472000Et02 434681E«01 412000E+02 


1 • 21 J » 11 (X.Y.2I . . .. 

472000E*02 . ..68614E+01 0. .72000E402 469527E+01 . 430146E+00 .... 472000E+02 470308E+01 _461140E*00 

472000E»02 . .70989E+01 493746E+00 472000E+02 .71586E+01 .. 412864E+01 472000Et02._. . 472110E+01 4l6641Et01 

.72000£*02 .72564E+01 ,20754F*01 .72000E+02 472946E*01 .25243E*01 . 472000E*02 473258E+01 430144Et01 

.72000E+02 473496E+01 ,, 43S475E+01 _ 472000E+02 . 473665E+01 441250E+01 472000E+02 473758E+01 447475E+01 

.72000E+02 .73772E+01 .54144E+01 472000E*02 , 473702E+01 . . 461245E*01 . . 472000E+02 473535E*01 468754E+01 - 

472006E+02 473250E+01 476641E*01 472000E*02 472807E + 01 484864E+01 . .72000E»02 .. . 472141Et01 _493375E*01 

472000E«02 .71137E^01 . .10211E«02 „ .472000E+02 .6939SE«01 4lll01E»02 . , 472000Et02 46714SE«01 4l2000Et02 


I - 21 J - 16 <X»Y»2) . . , . . 

472000E + 02 488751E+01 0, . „ ,72000E*02 „ _489577E*01 _.,430146E+00... 472000E+02 . 490334E+01 461140E+00. 

.72000E*02 491040Ei01 493746E*00.. 472000E+02 491703Et01 __4l2864E + 01 .72000E+02 .92323E+01 ,16641Ef01 

.72000£*02 .92897E+01 .. .20754E+Q1 .72000E402 493416E*0l 425245E+01 472000E+02 ..493872E*0l 430144E401 

.72000E+02 494253E+01 . .35475E401 .72000Ef02 494550E+01 ,41250E*01 »72000E*02_ «94750E*01 447473E+01 

.72000E+02 .94842E+01 ,34l44E+.0l., .._.72000E*02 494811Et01 46124SE401 .72000E+02 .94638E+01 .68754E+01 

.72O00E+02 .94291E+01 . ,7664l£*0l 472000E+02 493730E + 01 484864Et01 472000£*02 492890E + 01 493373E+01 _i_ 

.72000E+02 , 491685E + 01 ... .10211E*02 .72000E402 489991E#01 4lll01E>02 472000E+02 487646E401 412000E+02 



I - 21 J - 21 (X>Y#Z) 



,„.72000E+02 

. •U67eE402. 

. . 0. 

..72Q0QE+02_ 

__ai727E+02. 



-.72C00E*O2 _ 

....11817E402. 

.93746EtOO 

.W2000E.+ Q2... 

.11B59E*02. 

^ 

.. .72000E+02 ... 

...11939E + 02.. 

.2075<E*0L 

.^72000E + 02_ 

.^119756 + 02 


„ .72000E + 02 

__ a2036E + 02 

.35475E+01 

.720006 + 02.. 

.12058E+02 



_.72000E + 02 . 

.12083E+02 

..,.i5414.4E + 01 

372000E+Q2.- 

.12O02E + O2 


, .72000E402 

.12047E+02 

.76641E+01 

.72000E+02 . 

.. ...120086+02 

— 

_,72000E+02 . 

. .11069E+O2. 

a02HE + 02 

..72000E+02.. 

_..11760E+02. 








1 •. 21 J.» . 

.26. (X#y«Z)_ 





...72000E+02 

asoooEfOz. 

0. 

.^7200QE+02_ 

.15QQQE+02-. 


.72000E+02 

.1600CEiC2 

.•93746E+00 

..720006 + 02.. 

. — ..15000E+02 


... •72000E4U2 .. 

. a5000Ei02 

..... .20754E + 01 

.72000E+02 

.15000E + 02 


.72000E+02 . 

. a5000E+02 

.. .35475E + 01 

.720006+02.- 

.150006+02. 


.72000E402 

. a5000E402 

.54144E + 01 .. 

.720006+02^ 

.15000E + 02. 



. .72000E+02 

.a5000E402 

.76641E+01 

..72000E + 02_ 

_...15000E+02_ 

— 

.72000E+02 . 

. a5000E402 

.102116+02 

.72000E + 02... 

a5000E + 02.. 


I ■ 26 J ■ 

1 (X#Y>Z) 

- 



— 


•90000E+02 

.545896+01 

0. 7 

.900006+02 

' .54586E+01 


•90000E402 

•54560E+01 

...93746E + 00 

.900006+02 

- .54535E+01 


.90000E+02 

•54448E+01 

.20754E + 01 .. 

.900006 + 02.- 

.54380E + 01 



•90000E+02 . 

•54177E+01 

. .35475E+01. .. . 

.900006 + 02.. 

.. ...540326 + 01 


•90000E402 

.53629E+01 

•54144E+01 

.900006+02 

.53361E+01 


.90000E + CI2 

•52665E+01 

.76641E+01 _ 

.900006+02 

_„..52230E+01 


•90000E402 

.51172E+01 

.102116+02 

.900006+02 

.. .50550E + 01 

-- 

I - 26 J • . 

6* (X>YiZ) 

, . 




•90000E402 

•61275E+01 

0. 

.900006 + 02 ~ 

7 ” .’61273E + 01 


.90000E402 

.61249E+01 

.937466+00 

.900006+02 . 

•61225E+01 


•960UOE402 

•61144E+01 

.207546+01 

.900006+02 

.61O01E + O1 


•9000UE402 

.60892E+01 

.354756+01 

.900006+02 . 

.607576+01 


•90000E*02 

.60383E+01 

.54l44Ef01 

.900006+02 

.60133E + 01 


•9O000E402 

.59486E+01 

.766416+01 

.900006 + 02 - 

_ .59082E+01 



>90000E402 

.56098E+01 

.102116 + 02 . 

.900006+02 

.. . .57519E + 01 


I • 26 J - 

.90000E + 02 

11 (X>Y»Z) 

.72877E+01 

0. 

.90000E + 02. 

.728756+01 

.9000CE+02 _ 

. .72854E + 01 

. ...937466 + 00. 

.90000E + 02_ 

.72834E+01- 

.90000E+02 

.72763E+01 

.207546+01 

.900006+02 . 

... .727096 + 01 

•90000E+U2 

.72544E+01 

.354756+01 

.900006 + 02 . 

.72427E + 01 

.90000E+02 

.72101E+01 . 

.541446+01 

.900006+02 

.71885E+01. 

.90000E + 02 

.71322E+01 

.766416+01 

.90000E + 02_ 

.709706 + 01.. 

.90000E + 02 

.70115E+0X 

.102116+02 

.900006 + 02.. 

.69612E + Q1 


•301A6E+0a .72000£f02 .11T73E + 02 .611^0E*OQ 

.12864EtQl *,720Q0EtO2 ^900E.+ 02._ •16641Ei03 

.25245Ef01 ^72000Et02 .12007E+02 ^301^^E+01 

.41250E*01 .72000Et02 ^1207^E>02 .47475Ei01 

.6l2^5Ef01 .72000Ei02 ^12071Et02 •6875^E+01 

. .8486^E + 01 .72000E + 02 .11950E + 02 ..93375E+01 

.11101Ef02 .7200OE+02 ai618E+02_ .12000E^02 


^30U6EiQQ .72DO0£tO2 .1500DEf02. ^61140E-*00 

-_^a286^£ + 01 .72000E+02 .. ,>..15000E + 02 .16641E+01 

aSZ^SE+Ol- •72000E+02 •15000E402 - --..30144E+01 

_.41250E + 01 .72000E+02 *15000E + 02--l-.*47475Ef01 

•61245E + 01_.-..72000E+02 .15000E + 02 __.68754Ef01 

. .. _.84864E+01 .72000E^02 .15000E + 02 •93375E+01 

.lllOlE + 02 aZOOOE+OZ •15000E402 •12000E*02 


•30146E^00 

a2864E^01 

•25245Ef01 

•41250E+01. 

•61245E+01 

•84864E^01 

•11101E^02 


•90000Ef02 

•90000Et02 

•90000E402. 

•90000E^02 

•90000E402 

•90000E^02 

•90000E^02. 


-•54577E+01 

-•54490E+O1 

•54291E+01 

•53851E401 

•53041Ef0l 

•51732E+01 

•49868Ef01 


•61140EfOO- 

a6641Ef01. 

•30144E^01. 

.47475Ef0l 

• 68754E+01.. 

•93375E401 

•12000Ef02 


•30146Ef00 

•12864E^01 

«25245E401 

•41250E^01 

•61245E401 

•848646^01 

•I1101Ei02 


•90000Ef02 

•90000Ef02 

•90000E402 

•90000E^02 

.90000Ef02 

•90000E^02 

•90000£f02. 


•61264Ef01 
•61191Ef01 
•60999E^01 
•60569E^01 
•S9636Et01 
• 5e619Ei>01 
.a6666Ef01 


•61140E400 
•16641E401 
•30144E401 . 
..47475E401. 
.68754E401- 
•93375E401 
a2000E402, 


.30146E400 .90000E402 •72867E + 01 .61140E400 ; 

a2864E401 •90000E402 •72804E + 01 a6641£401 

• 25245E401 . aOOOOE + OE .72637E401 - •30144E401 , 

-.41250Ef01- .90000Et02 •72281Et01 •47475E+01 

- •61245E401,..-- .90000E + 02 . - .71626E401 .68754E401 I 

.84864Ef01 .90000E402 .^.70568£ + 01 - . .93375E401 I 

ailOlE + 02 .90000E402 .69062E401 ♦12000E+02 • 


,90000E^02 , 
.90000Et02 
•9000UE+C2 . 
.90000E’f02.... 
•90000E+02 
•900U0E+02 . 
..90000E + Q2„„ 


. ■ 26 J .• . 

.90000E+U2 
.90000E+02-. 
.90000E+02 . 
•90000E+02 . 
.90000E + 02 .. 
.9GOOOE+02 
•90000E+02 


■ 26 J ■ 

.90000E+02 
•9U000£^02 
.90000E+02 
•90000EiC2 
•90000E^02 
•90000E+02 
•90000E^02 


„ . .91960E + 01 . .0. .90Qa0Et02 a9195aEi0l ^3QlA6E*00_-_ ^9000QE♦0^ .9I952Et01 

.9l9't2Ef01 .937^6E>0Q ..9000QEt02 ^919^7Eta3 .I2B64Et01 .9OQ00E*O2 .9190^£tQl 

„....91874e + 01 .2075^E*Q1 ^90000E*02 .91033E+O1„ __*252.^5E+OX .90000Ef02 .91779E^0l - 

.91709E + 01 •35475E±OJ ..90000E+02 •9l62lEiOJ *.41250E>01 .900006*02 ^91511E*01 

. .913766*01 . . .5^1<r^Ef01 .900006*02 .912126*01 .... .612456*01 .900006*02. ..910186*01 .... 

__.907b9E*01 .766^16*01 .900006*02 .905246*01 .848646*01 .900006*02 .902226*01 

.890816*01 .102116*02 .900006*02 .E9502E*0r . 111016*02 .900006*02 .890886*01 


-.21. (XjIjZ) - 

.-.118526 + 02.- 0. .90000E + 02 .11852E + 02 .30146E + 00 .90000E+02 .118526+02 

._...11851E + 02 .93746E*0Q__^90000E*02 .118506*02 .120646*01 .900006*02 .118496*02 

. •116466 + 02 .2075^6*01 .900006*02 .118456*02 .252456*01 .900006*02 .118426*02 

__.11839E*02 .354756*01 .900006*02 .118346*02 .412506*01 .900006*02 .118286*02 

.11821E + 02 .. .541A4E + Q1 _.90000E*02 .118126*02 .612456*01 — .900006*02 .118016*02 

.117896*02 . - .766UE+01 .900006*02 .117746*02 .848646*01^^^ .900006*02 .117586*02 

.117396*02. .102116*02 .900006*02 *117196*02 .111016*02 —-.900006*02 .116966*02 


26 (X#Y#Z) . 

.150006*02 0. .900006*02 .150006+02 .301466*00 .900C0E + 02 .150006*02 

.150006*02 .937466*00 .900006*02 .150006*02 .-.128646*01 .900006*02 ^150006*02 

.150006*02 .207546*01 •900006*02 .150006*02 .252456*01 1 . .900006*02- .150006*02-^. 

.150006*02 .354756*01 . .900006*02 .150006*02 .412506*01.. .900006*02 .150006*02 __ 

.1500CE + 02 .541446*01 .900006*02 .150006*02 .612456*01 .900006*02 .150006+02 

.150006*02 .766416*01 .900006*02 .150006*02 .846646*01 .900006*02 .. .150006*02 

.150006*02 . .102116*02 .900006*02 .150006*02 .111016*02 .900006*02 .150006 + 02 


-.611406*00 - 
-.166416*01- 
-.301446*01 
..474756*01 
•687546*01- 
_.93375E*01_ 
-120006*02 


.611406*00- 
-.166416*01 
^301446*01- 
.474756*01— 
.687546*01 
.933756*01 
.120006*02. 


.611406*00 
-.166416*01 
•301446*01 
-.474756*01 
.-.687546*01. 
• 933756*01. 
•120006*02 


III. Conclusions and Recommendations 

A method for numerically generating three-dimensional 
coordinate system for a wedge with a curved after body has 
been presented. The computer program AFTBDY used to numerically 
generate the coordinate system stores the data on a disk file. 

This coordinate system will become a part of the input to the 
three-dimensional Navier-Stokes computer program CANS (developed 
by Dr. Julius E. Harris at the NASA Langley Research Center) . 

The computer program CANS is a code to solve three dimensional, 
time split, viscous, compressible Navier-Stokes equations. 

No Navier-Stokes computer runs using CANS have been made as 
yet for a wedge curved after body using the coordinate system 
generated by the program AFTBDY. There are tv/o areas in program 
AFTBDY where improvements may have to be made in case the Navier- 
Stokes solution generates errors which are attributable to the 
coordinate system. 

The first area could be the way the attraction term Q is 
computed. Q is presently computed based on the y variation on 
5=1, 5=1, n = 1, JMAX line. For all n planes the same value 
of Q is used. However, as the total y variation from ri = 1 to 
n = JMAX for various 5 and ? stations changes, this treatment of 
Q gives ri coordinate line variation as shown in Figure 10. Close 
to the after body, the n lines do not follow the after body contour 
satisfactorily. This problem can be minimized by choosing Q as a 
sum of exponential functions as suggested by Thompson, Mastin, and 
Thames of the Mississippi State University. 
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The second problem area is the use of Dirchlet boundary 
conditions on ? = 1 and ^ = KMAX planes . This results in abrupt 
change of slope near the ^ = 1 and C = Kr4AX planes (Figure 10) . 
This problem can be avoided by using Neumann boundary conditions 
on these two planes instead of Dirchlet boundary conditions . 



IV, Proqrammer/Analyst Section 


Program AFTBDY is structured such that the array size in 
^,r\, and 5 direction can be varied from problem to problem by 
the use of text editor on the CDC 6600 machine. The array sizes 
are coded to be IMXLL, JMXLL, KMXLL in rif C direction, respectively. 
These are changed to the appropriate numeric number by the text 
editor prior to compiling the program. The job control cards, input 
explanation and other helpful hints are provided towards the beginning 
of main program AFTBDY. The program is extensively commented. 

Though, the code is not structured in the true sense, extreme care 
was taken to avoid unnecessary upward branching in the code. A 
brief discussion of each subroutine and their function follows . 


A. Main Program AFTBDY 

This program performs the following function, 

i) Initialize variables to zero. 

ii) Read in input data. If the field size is larger than the 
maximum field size allowed then the program aborts with a 
diagnostic message. The maximum size allowed is ILIMIT, 
JLIMIT, KLIMIT and field size input is II/IAX, JMAX, KIIAX. 
ILIMIT, JLIMIT, KLIMIT are set by the text editor. The 
storage in N.S. code is such that it is necessary to set, 

ILIMIT = IMAX 
JLIMIT = JMAX 
KLIMIT = KMAX 

iii) Call subroutine INIT to specify initial guess and compute 
P, Q, and R. 

iv) Call IPRTC to print initial guess, 

v) Call IPLTC to plot initial guess. 
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vi) Call CALCOR to iterate for the value of y CALCOR passes 
back flag KKK = 1 if the solution converges, otherwise 
KKK remains set to 0. 

vii) Call IPRTC to print partially or converged solution. 

viii) Call IPLTC to plot partially or converged solution. 

ix) If the solution converges , then store solution on unit 
1. If the convergence is not reached in ITERI4X itera- 
tions, then a diagnostic message is printed and the 
solution is not stored on unit 1. 

B. Subroutine INIT 

This subroutine computes the initial guess and the inhomogeneous 
terms P, Q, and R. 

C. Subroutine CALCOR 

This subroutine computes for y. A do loop is set to perform 
ITERMX iterations. Maximum error in y and its I, J, K location 
is stored in YERRMX, lYERl, lYERJ and lYERK. After each iteration 
maximtim error is checked against the maximum error allowed. If 
maximum error is less than maximum error allowed, then the flag KKK 
is set to 1 and an exit from the do loop is taken. If convergence 
is not reached in ITERMX iteration, then the flag KKK remains set 
to zero. 

D. Subroutine IPRTC 

This subroutine prints initial guess, partially converged 

solution or the converged solution. The following flags indicate 

the status of the solution 

ISOLN = 0 Initial guess 

ISOLN = 0 and KKK = 0 Partially converged solution 

ISOLN = 0 and KKK = 1 Converged solution 

In case of ISOLN =0, IPRTl flag indicates whether the 

solution is to be printed or not. IDELl and JDELl indicate the 

interval of ? and n at which the solution is to be printed. 



When ISOLN =1, IPRT2 flag indicates whether the solution is 
to be printed or not. In this case IDEL2 and JDEL2 take on the 
same significance as IDELl and JDELl in the case of the initial 
guess. 

E. Subroutine IPLTC 

This subroutine is used to plot initial guess, partially 
converged solution or the converged solution. The flags to 
indicate the type of solution are the same as in the subroutine 
IPRTC. This subroutine plots 5 and C planes as indicated by 
input to the program. 

F. Block Diagram 

The block diagram is presented on the next page. 
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F. Block Diagram 






TO PLOT 
INITIAL GUE 3S 
(ISOLN = 0) 



:all iprtc to 

'RINT PARTIi iLLY 

lONVERGED ( ISOLN = 1, KKK =0) 

)R CONVERGE ) (ISOLN = 1, KKK = 1) SOLUTION 



61 





21 


AFTER BODY NUMBER 1 
TEST STARTED FEB 16,1979 

— INITIAL GUISS — 

IMAX 26 JMAX 26 £MAX 

- y-Z PLOT I 1X0. 





INFOPLT 1 


Figure 3. Initial Guess 5=1 Plane 









AFTER BODY NUMBER 1 
TEST STARTED EEB 16,1979 

— INITIAL OUISS — 

IMAX 26 JMAX 26 KMAX 21 


X-y PLOT X 


0 . 


Figure 6. Initial Guess ? = 1 Plane 
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AFTER BODY NUMBER 1 
TEST STARTED EEB 16,1979 

— INITIAL GUSSS — 

IMAX 29 JMAX 26 KMAX 

- X-Y PLOT X 11 Z .41260X 01 





Figure 7. Initial Guess ^ = 11 Plane 



AFTER BODY NUMBER 1 
TEST STARTED FEB 16,1979 
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Figure 8. Initial Guess ? = 21 Plane 
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APTER BODY NUMBER 1 
TEST STARTED FEB 16,1979 

— aONVXa.GlD SOLUTION — 

IHAX 26 JMAX 26 XMAX 21 

~ Y-Z PLOT I 1X0. 


aJO 




INFOPLT 7. 

Figure 9. Converged Solution 5=1 Plane 
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AFTER BODY NUMBER 1 

TEST STARTED FEB 16,1979 

— aoNyaaoxD solution — 

IMAX 26 JMAX 26 XliAX 
-- y-Z PLOT 1 22 X .766001 02 



X 


3 > 

INFOPLT 8. ^ 


Figure 10. Converged Solution C = 22 
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AFTER BODY NUMBER 1 
TEST STARTED EEB 16,1979 

— GONYEB.GXD SOLUTION — 

IMAX 26 JMAX 26 XMAX 

» Y-Z ?LOT 1 26 X .90000E 02 
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AFTER BODY NUMBER 1 
TEST STARTED EEB 16,1979 

— GONVIRGID SOLUTION — 

IMAX 26 JMAX 26 KliAX 21 

-X-yPLOT X 1 Z 0, 




Figure 12. Converged Solution C = 1 Plane 
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AFTER BODY NUMBER 1 
TEST STARTED FEB 16,1979 

— QONYKilClD SOLUTION — 

IMAX 20 JMAX 20 KMAX 21 


- X-Y PLOT S XI Z .412601 01 




Figure 13. Converged Solution c, 


11 Plane 



ABTER BODY NUMBER 1 
TEST STARTED FEB 16,1979 

— aONYSHGlD SOLUTION — 

IMAX 25 JliAX 26 XJiAX 21 

- X>Y ?LOT X 21 Z .12000X 02 





Figure 14. Converged Solution ? = 21 Plane 
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